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One of the great rewards enjoyed by persons who spend their 
lives in scientific research work is the unpredictability of 
what they are doing. It is an exciting world and one I enjoy 
immensely. However, not all of the surprises are the kind one 
likes to boast about. Principal Investigators have to admit to 
less productive years as well as proudly present their successes. 

On this program, in the meteorological study area, we have had 
one of those years. While good work has been done by several 
persons, it is as yet too incomplete to include in this report. 

We expect to present these efforts in next year's report, 

I am happy to present three papers by Dr. Aniruddha Das 
and his principal advisor. Professor T. C. Huang. Publication 
of these papers concludes Dr. Das' development of a generalized 
flexible satellite attitude control model and the application 
of that model to some relatively simple analyses. We anticipate 
that Das 1 model will be used by government agencies and by 
industry in more complex applications. 

I am especially grateful to Professor Huang for his assistance 
and support. We sincerely appreciate the patience and support of 
the many dedicated persons in the National Aeronautics and Space 
Administration with whom we have worked during the past year. 

Verner E. Suomi 
Principal investigator 
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STABILITY OF STOCHASTIC SATELLITES 


T. C. Huang and Aniruddha Das 


ABSTRACT 


The effects of random environmental torques and noises in the moments of 
inertia of spinning and three-axes stabilized satellites are compared analyti- 
cally and by analog simulations. Four analytical methods are used to compute 
the mean values and variances of the satellite response. Among the analytical 
methods, it is shown that the Fokker^Planck formulation yields predictions 
which most coincide with the simulation results. The variances of the responses 
have been shown to have an initial period of growth. This growth rate falls 
off with time and the variances reach and stay at an equilibrium value. The 
growth rate is also shown to be an increasing function of the inertia noises 
and the nominal spin rate. 


NOMENCLATURE 


A ± , i = 1-4 
a i , i - 1-27 
C 

D l* D 2 
F.{F 1 } 

f*,f*(w,t|m,T) 


JUJfjh i - 1,2,3 
I, (Ip 5 i - 1,2,3 


I 1, I 2 , I 3 

T r T 2’ T 3 


K 

L,L(6i, 02,6^|u,t) 


■ Arbitrary constants; Eq. (74). 

* Coefficients defined by Eqs. (10-18) and Eqs. (19-27) . 

* Arbitrary constant; Eq. (74). 

* Arbitrary constants; Eq. (74). 

- Vector forcing function; Eqs. (91, 92), 

■ Conditional joint probability density function of m(t) 
given the values Of j3(t) . 

* Arbitrary random forcing functions; Eqs. (1) , (19)-(2l) . 

* Mean values of f ,{f^}. 

* Components of M 20& , M 0 20* ^110* res P ectivel y 5 Eqs* (74), 
(74a) , (74b) and (74c). 

* Stochastic moments of inertia of the satellite; Eq. (1) . 
■Mean values of 1^, I 2 and X^, respectively, 

■Functional defined by Eq. (95) . 

* Polynomial function of p ; Eq. (70). 

■ Derivative characteristic function with parameters 8^, 

Q 2 and .8 j for the random variables for a given c3(t) ; 
Eq. (7). 


1 


2 


L* 

H ij5 i,j - 1-6 
“kirn 

N 1 , ± - 1-7 
N ij; i,J - 1-6 

r 

r ±j J i,j - 1-4 
* ± y i.j “ 1-4 
T 

t 

V i.J - 1-* 

u, { Ui }; i - 1-6 

v, {v ± }; i - 1-6 

V 1 - 1-8 

a iJ ’ i ■ °- 6 
e ± , i - i-3 
By, j - 0,1, 2,3, 
«(t) 

« ± , i - 1-3 

e 

e i t i - 1-3 
D i , 1 * 1-3 
e i , t - 1-3 
Xf ±> i « 1-3 

Af t , 1 « 1-3 


Matrix differential operator; Eq. (79). 

Covariance matrix of u; Eq. (5). 

Statistical moments of w(t) for a given w(0) ; Eq. (30) . 

Parameters related no by Eq. (108). 

Covariance matrix o': v; Eq . (94) . 

Eigenvalue of various equations. 

A measure of the noise levels; Eq. (122). 

Coefficients defined by Eqs. (74e) - (74g) . 

Coefficients defined by Eqs. (74e) - (74g). 

Period of time in which the most-likelihood estimates 
of are required. 

Time. 

Coefficients defined by Eqs. (74e) - (74g). 

Random vector; Eq. (4). 

Random vector; Eq. (93) . 

Coefficients of the characteristic polynomial for p; 

Eq. (41). 

Components of a^; Eqs. (46), (58) etc. 

Lagrangian multipliers; Eq. (95). 

Components of 6^; Eq. (109). 

Dirac's delta function. 

White noises associated with A'f^; Eq. (2). 

Largest absolute value of f° r a H 1 ant * j» (108a) . 
Sample space white noises associated with 1^; Eq. (75). 
Time dependent white noises associated with 1^; Eq. (75). 
Parameters of L; Eq. (7). 

Total forcing functions defined by Eqs. (10) - (12) and 
Eqs. (19) - (21). 

Mean values of Xf^. 
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X'f lt i - 1-3 

■ Total forcing functions defined by Eq. 

a). 

X’f t . i - 1-3 

* Kean values of X’f^. 


X 1 ,X 2 

* Parameters defined by Eqs. (71), (72). 


V ± , i - 1-3 

* Total white noises associated with 1^; 

Sq. (2). 

*1 

* Parameter defined by Eq. (74d) . 


P 

* Parameter defined by Eq. (74d). 


p kto 

“ Statistical coefficients defined by Eqs. (6), (8) 

° v i - 1-3 

= Standard deviations of Eqs. (115), 

(116). 


* Nominal spin rate of the satellite. 



* Nominal angular velocity vector of the 

satellite. 

<*), 

* Angular velocity vector of the satellite; Eq. (1) 

w 9 l» ± } 

■ Realized angular velocity vector corresponding to 

j * 1-34 

* Components of Eq. (76). 


OPERATORS 



E{ } 

* Statistical expectation. 


o 

* Mean value. 


[ ] T 

= Transpose. 


o 

d 

* = dt ’ 


INTRODUCTION 




This study compares the effects of stochastic geometry and random environ- 
mental torques on the pointing accuracy of spinning and three-axes stabilized 
satellites. A comparison of pointing accuracies requires a comparison of the 
rates of error growth over and above any criterion for the asymptotic stability 
of the satellites. For this reason, this study is oriented towards the deter- 
mination of the statistical properties of the Satellites ’ responses. The 
questions of stability have been answered indirectly by the computed responses. 

The reason for considering the environmental torques on the satellites 
as random is self-evident. The geometries of the satellites are considered 
stochastic in order to have a phenomenological model of the motions of the 
satellites’ flexible structural elements. If a Satellite were absolutely rig- 
id , its inertia properties would have been constant for all time and measured 
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to a near certainty. Because real satellites contain many flexible and moving 
parts, their moments of Inertia can be assumed to be stochastic variables vlth 
certain associated noise. 

To be more specific, the rigid body Euler's equations 
hh + (1 3- 1 2> w 2“3 “ X ' f l 

1 2“2 + (I r 1 3 ) “l“3 “ X ' f 2 (1) 

Vs + (I 2“ I l ) “l“2 ’ X ' f 3 

governing the motion of satellites will now be analyzed. In the above equation, 

I- , I„, I~ are the stochastic principal moments *of inertia of the satellite. 

± A 3 T T 

The vectors jd * [ (0 l* u 2 >tt 3^ and “ [X'f^rA'fjjA'fj] are the angular velocity 

vector and the environmental torque vector of the satellite, respectively, along 
the principal axes of inertia. And X' is a parameter. The vector X ? JE and, 
consequently, the vector jd are random variables. 

Equation (1) is an example of an intrinsically nonlinear system of equa- 
tions with random coefficients. The difficulty of obtaining an explicit solu- 
tion to Eq. (1) can be appreciated when we realize that the stochastic version 
of even a simple scalar linear equation is actually nonlinear due to the de- 
pendence of the solution on the random coefficients. (See Refs. 1, 2.) The 
situation has been made even more complex by the presence of several contradic- 
tory methods for solving stochastic equations [1], A widely used method of 
solving stochastic equations is the Fokker-Planck approach. In this, the equa- 
tions are assumed to define a Markoff process and the transition probability 
densities of the responses are computed directly as a function of time. Sev- 
eral interesting equations have been solved by this method in Refs. [3-7]. 

Another useful method, using perturbation techniques for solving stochas- 
tic equations, was discussed in Refs. [8,9]. This is one of the "honest*' meth- 
ods in which response is solved analytically in terms of small random param- 
eters. The stochastic properties of the response are obtained from the analy- 
tic solution as secondary results. 

A third promising method of solution can be obtained by extending the line 
of logic shown in Ref. [10]. This method determines the most likelihood estimates 
of the response by maximizing the joint probability density of all the stochas- 
tic variables of the system. This is essentially a formulation of the Kalman 
filter for the case of deterministic coefficients and random forcing functions. 

Lastly, there is the obvious method of initially assuming the system of 
equations to be deterministic and then attributing the proper stochastic prop- 
erties to the deterministic solutions. It is, of course, true that this method 
is rigorous only If the random parameters are constants in time. The stochastic 
properties of the eigenvalues and eigenvectors of such systems have been comput- 
ed in Refs. [11,12]. This method is worth investigating for slowly varying 
parameters with random step increments. 
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The response vector, a>(t), of a rigid satellite governed by Eq. (1) 
will be analyzed using the above mentioned techniques. The analytical responses 
are then compared with results of an analog computer simulation. This allows 
verification of the relative merits of the analytic methods. 

THE FOKKER-PLANCK APPROACH 

This method of obtaining the response characteristics of stochastic equa- 
tions is based on the analysis shown in Refs. [1,13] . The application of this 
method on Eq. (1) proceeds as follows: 

Let the random variables y 2 , y 3> 6 ^, $ 2 and 6 ^ be defined by the equa- 
tions 

I - X + y 8 1 - 1,2,3 

( 2 ) 

X*f -7 f I ± + 6 ; 1 *1,2,3 

The bar on top of a symbol indicates mean values. Hence, 

*M± -7 ± « ° 5 1 - 1,2,3 (3) 

Let the stochastic vector u be defined as 

U- [y 1 ,v 2 ,P 3 ,X’f 1 ,X , f 2 ,X , f 3 ] T (4) 

It is assumed that y^ and 6 ^, i * 1-3, are white noise disturbances, such that 
the matrix elements, , i,j = 1 - 6 , are defined by 

ECujU } * (5) 

In Eq. (5) and in the following, (t) is the Dirac’s delta function and the 
operator E{.} denotes statistical expectation. 

Let P^ m ( a3 i» U) 2 ,w 3 , t) be the statistical coefficients of various orders 
where are the realizations of the responses for i * 1-3, at any point 

in the time and sample spaces. Let It also be defined that f*[w,t |w(0) ,0] 
is the joint conditional probability density of the response vector, w, given 
the values of <S(0) at t * 0. Thus, 

P U m - At*) 1 1 1 <V“l )k( “2"“2 )4(w 3 _ “3 ),l,£ * < ^* t+At ^' t)d “l d “2 d “3' (6) 


Although Eq. ( 6 ) is used to define the coefficients these are usual- 

ly calculated from the derivative characteristic function 


L(8^»9 2 , 63 j t) 


This, in turn, is defined by 
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Ue v 6 2 ,6j\!£.,t) - jjjjj ^ [E{exp[i Sj { Uj (t+At)- Uj (t) }] |w,t} - 11 (7) 

vhere i ■ /-I . 


Comparing Eqs. (6) and (7), an alternative definition of p # can be 
obtained as 


. -(k+i-m) . a* 1 *** l . 
kto 3e k 30*96^ 


6 ! ‘ ®2 ■ fl 3 ■ 0 


( 8 ) 


Let It be assumed that 


ojj (t+At) - Uj(t) * Wj(t)-At 


(9) 


The values of are now easily calculated from Eqs. (7,8,9). For example, 


. r 3L, 
P 100 ” " ll 36 1 J 


e, - o * - 1 AtHo [E{exptl A e j“j Atl 


e 2 - o 
e 3 -o 


le.-o 
le 2 - o 
|S 3 - 0 


or P 100 “ at-*) TE [E{ “i At l“.t>] 


or P 100 ‘ fit^O It [E{ T^T [X ’ f l “ (I 3 - I 2 +p 3 - p 2 ) “ 2 “ 3 1}3 


; at^O At E{ f (1 ' l' )lX ' f l ‘ < l 3 -I 2 + ‘ , 3 _u 2>“2“3 ]} 

1 L 1 

Expanding the right hand side and neglecting the cubic and higher order terms 


in y 


r 


M n -Mi, A A M 1A __ 

[ ( % 12 + I 0 - T 3 }a> 2 o ) 3 - + X’f , ] 


K ioo I x 

1 1 


Proceeding similarly, it is easily seen that 


( 10 ) 


p 010 ■ f l{ ~^ + *3 - * 1 * 1*3 - # +X f 2l 


(ID 


y 001 


1 _ _ A A 


( 12 ) 


7 


4(1 -f,) 

p 200 “ y2 ^ {M 33 " ^23 + M 22 " 17 (M 13 -M 12 )1 “2“3 


+ {j- (I 3 -I 2 )M 14 - 2(M 34 -M 24 )}^ 3 + M 44 ] 


(13) 


p 110 • [{M 13 + M 23 - M 33 - M 12 + Vli) 

(I, -I,) (I,-I 2 ) do -I, ) 

+ -l-d- (M 13 -M ) .+ (M 23 -M 12 ) + -4-d- (M 23 -M 22 ) 

I x A 2 *2 

(I -T )(l -i ) . o d 3 -i ) 

“ • -- T^ ' M 12 } “l“2“3 + {M 34 - M 14 ” M 14 


(1,-1, ) (1,-1,) 

^=d- m 24 }£ 1 £ 3 + {m 25 - m 35 + - A =— m 15 


<W 


T M 25 } “2“3 + M 45 ] 
2 


(14) 


(Xo-I 0 ) 


p 10l“iY t{M 13 “ M 23 ” M 12 + M 22 + \ (M 11 _M 12 ) 

I 1 I 3 1 


(1,-1,) (1,-1,) (i,-x,) 

“V“ VlZ 5 + ~T (M X3‘ M 23> “ ~^T (M 33- M 23^ 

L 1 3 3 


<vv<w 


„ 2 ~ 


< I 1- I 2> 


f Y ~ ~ M 13 } “l“2“3 + {M 14 - M 24 - “ Y- M 14 

13 1 

(1,-1,) (T,-I,) 

- -V" M 34 }£ 1“2 + {M 36 - M 26 + J — M 16 
3 1 

( 1 ,- 1 ,) 

+ -^-« 3 6>“ 2 “3 + M 4 6] 

3 


( 15 ) 


8 


“020 ’ 3 t {M 33 " *13 + M 11 - f (I 3" 1 1 )(M 23- M 12 )} “l“3 
X 2 2 

+ {2(M 35 -M 15 ) -f- (x 3 -i 1 )M 25 }£ 1 £ 3 + M 55 l 
2 


(16) 


“oil “ II 1 -13 " M 23 " M 11 + M 12 + I 1 (M 22" M 12 ) 
2*3 i 2 

(1,-X,) (1,-1,) (1,-1,) 

+ ~J~ (M 23- M 12> + (M 23- M 13 ) + T; 

a 2 3 3 


(1,-1, ) (1,-1,) , (1,-1,) 

- TI M 23 } V2“3 + <«36 ' M 16 * M 26 


T T a l i 3 jb ib i 

2 3 2 


- — 3 - — M 36 } “l“3 + {M 15 “ M 25 + 2 x X M 25 
*3 2 


<V I 1> 


M 35 } “l u 2 + M 56^ 


(17) 


M 002 


> i ( * 2 r, 2 


_ 2 [(M 22 - 2M 12 + M i;l - (I 2 -I 1 )(M 23 -li 13 ))ii) 1 0) 2 


+ { f- (x^W* - 2 (m 26 -m )}£ £ 2 + m 66 ] 

3 


(18) 


All of the first and second order expressions of are listed in Eqs. (10- 

18) above. The third and higher order are usually small and can be neglect- 

ed . Suitably defining the set of constants a , J « 1-27, Eqs. (10-18) can be 
rewritten as J 


“l00 “ a l“2“3 + Xf l - a 2 


p 010 ” a 3“l“3 + Xf 2 - a 4 


p 001 “ a 5“l“2 + Xf 3 " a 6 
p 200 “ a 7“2“3 + a 8 5 2“3 + a 9 


(19) 

( 20 ) 
( 21 ) 


( 22 ) 
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p 110 ■ a 10“l“2 S 3 + a ll“l“3 + a 12“2 S 3 + a 13 (23) 

p 101 - a 14VA + a l5V2 + a X6“2 S 3 + a 17 < 24 > 

p 020 - a 18“l“3 + a 19“l £ 3 + a 20 (25 > 

p 011 ’ a 2lW3 + a 22V2 + a 23“l“3 + a 24 (26) 

p 002 ' a 25“l“2 + a 26“l £ 2 + a 27 (27) 

Because the values of , corresponding to the system given by Eq. (1) 

are at hand, the Fokker-Planck equation involving the density f*[w^t | 5(0) *0] 
for that system can now be set up. This equation for the density is [I] 

3f* y .^vk+Ji+m k+£+m 

3t k+£+m>0 ( kU!ra! TXTTIm (p kAm f * } * (28) 

Substituting Eqs. (19-27) in Eq. (28) and neglecting all third and higher order- 
derivatives, Eq. (28) reduces to 

3£* 1 , -2-2 . - - . . 3 2 f* . 1 ' -2-2 ' a a . ,9 2 f* 

It" “ 2 Ia 7“2“3 + a 8 u, 2“3 + a 9 ] + 2 ta 18 u l“ 3 + a 19“l“3 + a 20 ] 3T" 

3a)^ 9 o>2 

2 

, 1 r -2-2 4 a * , t 3 f* , r ft a ^2 , - 

+ 2 la 25“l“2 + a 26“l u, 2 + a 27 ] + [a loV2“3 + a ll“l“3 

30)3 

2 

+ a 12 S 2 £ 3 + a 13 ] -j£^ + ivA + a 15 “l“ 2 + a 16“2 i 3 

. 3 2 f* . r -2- - . - - . - - . , 3 2 f* 

+ a 17 ] 33^ + Ia 2l“l“2“3 + a 22“l“2 + a 23“l U 3 + a 24 ] 3S^ 


2 2 — — 3f * 

+ [a 10 S l“3 + a 12“3 + a 14“l“2 + a 16 fi 2 ' W3 * X£ 1 + a 2 J 3^ 

7 2 A ^ . 3f^t 

+ [a 10 fi 2 2S 3 + a xl u 3 + + a^ - ayu^ - 

+ [a u 3 2 £ 3 + a 15 S 2 + 32^3 + a^ - a^^ - Xf 3 + a fi ] 

+ ta 10 £ 2 + a 14 S 2 + a 21 5 2 ]f* . 


(29) 
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The values of the density function can be obtained by solving this for- 
midable linear second order partial differential equation. But little useful 
information is obtained from the density function. The truly useful statisti- 
cal parameters are the mean values, variances, covariances, and other higher 
order moments of the satellite response* These parameters form a family, ML _ , 
which is defined by * 


aIc a % g A 


* / / . / f*[u l t;|u(0) r 0]du 1 d(a 2 dm 3 


(30) 


_00 —00 _oo 


and hence 


« / / /i^2 ( V n ) m d “ 1 d‘ 5 2 d “3 • <31) 


where R is the nominal value of the spin rate. Substituting the expression 
d€ * 

for from Eq. (29) into Eq. (31) and integrating, it is seen that 


M 100 ■ a l°*0l0 + a l M 011 - a 2 + Xf l 


(32) 

(33) 

(34) 

+ a 8*011 + a 9 (35) 
*110 “ (a ll° + X *2 " a 4 )d 100 + (a 12 !M '^^r a 2 ) *010 + a l°^020 + a 10 n *110 


*010 ” a 3 R *100 + a 3*101 " a 4 + Xf 2 
*001 “ a 5*110 “ a 6 + X *3 

*200 “ 2(X ^l -a 2 ) *100 + a 8 n *010 + 2a l JI *110 + a 7^ 2 *020 


+ a ll*101 + a 3^*200 + a 12*011 + a 13 
*101 ■ < X V a 6>«100 + a 16 n *010 + ( W 001 + a 15*110 


+ (a 1 SHa 16 )M 0U + a 17 


2.-. 


*020 “ a 19 tl *100 + 2 ^ Xf 2 -a 4^*010 + a 18° *200 + 2a 3 tl *110 
+ a 19*101 + a 20 

*011 ” a 23®100 + * X *3 -a 6**010 + (X ^2 _a 4 ) *001 + a 22*110 
+ ( a 23 +a 3 !J) ^ioi + a 2 4 
*002 ” 2(X *W*001 + a 26*110 + a 27 


(36) 

(37) 

(38) 

(39) 

(40) 


In deriving Eqs. (32-40), all third and higher order moments have been neglect- 
ed . Solving these nine first order ordinary differential equations, the mean 
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values, the variances, and the covariances of the satellite response are ob- 
tained completely. 


THE FOKKER-PLANCK RESPONSE 


At this point, it will be interesting to analyze the response predicted 
by Eqs. (32-40). These predictions will later be compared with an analog simu- 
lation of Eq, (1). 


^Let it be assumed that, at t *= 0, all second order moments (k+Jl+m - 2) 
and are equal to zero. In this stage, the satellite will behave as it 

does in the deterministic situation, that is, it will begin^to p recess with 
a rate proportional to ft. Then, as the values of M QOi and M q q 2 grow with time, 

the precessing rate and the nutation angle will also grow. Finally, the satellite 
topples down. This phenomenon occurs physically and in simulations. Thus, Eqs. 
(32-40) predict that the satellite response is greatly sensitive to the values 
of a^, (Af^-a^), a 2 g and a^. Because a^, and a^ are non-negative, these 

equations predict that an uncontrolled satellite governed by Eq. (1) is inher- 
ently unstable in the presence of random errors. The same conclusion can be 
drawn by applying the stability criteria of Refs. [14,15] to Eq. (1). The 
error growth rate of the satellite response can be minimized by minimizing the 
values of a^, a^, and Xf^. This can be done if * 0, = T 2 and the 

matrix is a diagonal matrix. 

The relative rates of error growth of spinning and non-spinning satellites 
will now be examined from the characteristics of the eigenvalues of Eqs. (32- 
40) . It can be shown that the eigenvalues of these equations satisfy a ninth 
degree algebraic equation of the form 

9 , 8 , 7 . 6 , 5 , 4, 3. 2, A ..... 

p + a Q p + a 7 p + ot 6 p + « 5 p + a^p 4- a 3 p + ct 2 p + ctjP » 0 (41) 

where a^, I s * 1-8, are appropriate constants. 


It is obvious that to have bounded growth rates, for all i must be non- 


negative. It can be shown that 
.2 


(X 3 -X 2 ) 


V= ^ - - TT lM 13 + M 23 - M 33 - M 12 + *IS*U ) 

A 1 A 2 A 1 


(1,-1,) (1,-1,) (1,-1,) 

+ (M 23 -M 12 ) + -4-i- (m 23 -m 22 ) 

1 2 l 2 


( 1 3 -I 1 )( I 3-I 2 ) 

V 2 


M 12 J 


( 42 ) 


Because usual satellite geometries are such that 
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j h >m * ll vh ] 


“ Eq. (42) says that cig ^ 0 if, and only if any of the following conditions exist: 

0 

1 

et 

(43) 

max[Hi 3 ,H 23 J < min ^ M 11 » M 1 2 * M 22 ’ M 33 ^ 

(44) 

In particular, ag >_ 0 if M ^ 2 >. 0 and 


M 13 - M 23 * 0 

(45) 

since min [M 11 ,M 2 2 *M 33 ] >_ 0 . 


Equation (45) states that one of the conditions 
rate is satisfied if the inertia noises in I- and t* 

for a bounded error growth 
are independent of the 

noise in I y But this condition usually is not satisfied because 

[ I 3 ■ ti + Ij and Ig ■ + I 2 and hence 


v 3 ~v 1 + v 2 


and , therefore , 


M 13- M 11 + M 12 


1 M 23 " M 22 + M X2 * 

I" 


l 

Thus, at this point it appears that Eq. (43) provides the only suitable con- 
straint and that this constraint is available only to three-axes stabilized 
l satellites. 

Now, let the conditions required to make a_ non- 
It can be shown that is of the form i 

'negative be considered. 

a 7 * °70 + a 7i a + a 72° 2 + a 74 flA > 

(46) 

where 


°70 " ^°22 a 12 + a 23 a 16 + *15*11^ 

(47) 

“71 “ ” 2la l a 23 + a 3 a 16 ] 

(48) 

“72 * - 6 a l a 3 

(49) 

“74 ' “ a 7 a 18 

(50) 

Another reasonable assuaptlon we can make now is that the inertia noises, V . , 
are independent of the forcing funcltons, Xf^. Assuming this. 
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a 2 " a 4 * a 6 " a 8 * a ll “ a 12 “ a l5 * a l6 


a 19 " a 22 * a 23 * a 26 * 0 


(51) 

(52) 


Using Eqs. (51,52), the criterion for non-negative ct^ becomes either Eq. (43) or 


“72 + a 74 B - 0 


(53) 


Equation (53) can be expanded to obtain 


6a l a 3 + a 7 a 18° - 0 


6 [ M 13" M 12 


hh h 


(I 3 -I 2 )][(I 3 -X 1 ) + 


M 12~ M 23, 


it 


A (1,-1,) 

_ 2 _ 2 [M 33 - 2M 23 + M 22 (M, ,-M, 0 ) ] [M. 

iti o i 


13 12' 


33 


2M, 


12 


12 


«<VV 

+ M U (« 23 - M 12 )1 


(54) 


Assuming the satellite geometry to be given by 
1 


I * I a T 
2 3 1 2 


V- Ui + v 2 


(55) 

(56) 


and that M are small compared to I,, Eq. (54) can be further simplified to 
read 1J 1 


2l l ~ M 11 (3M 22 ' 2M ll ) “ 2 - 0 


(57) 


Equation (57) is almost certainly satisfied for all real satellites and hence, 
cu is almost certainly positive. Equation (57) also states the obvious fact 
that , in the presence of inertia noise, a high spin rate tends to make the sat- 
ellite unstable. 


The expressions for will now be considered. It can be shown that ot^ 
is given by 


“6 “ “60 + “61 + “62 n2 + “63^ * “64^ + “66^ 


(58) 


where 
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°60 “ " ta X2 a X5 a 23 + a ll a X6 a 22 + ^*l" a 2 ) ^“XS^S®!!* 

+ <Xf 2 -a 4 ) ( a i a 22 +a 5 a l2> + ^^3 _a 6 ^ a 3 a 16” a X a 23 ^ ' 

“6X “ - [a 3 a 8 a 22 + a X a X5 a X9 + 2a 3 a X2 a X5 + 2a X a XX a 22 + 4a l a 3 (l£ V a ^ 

°62 “ a X0 a X6 a 23 
°63 ‘ 2ta X a XO a 23 + a 3 a XO a X6 J 
“64 “ 2la X a 3 a X0 " a 3 a 7 " VxS 1 
a 66 * a 7 a 10 a 18 

It has already been mentioned that, if (Xf^-a^) is non-zero, then even the de- 
terministic response is unbounded. Hence, to make any useful comparison, it 
must be assumed that flf^-a^) zero or ha» been made so by appropriate 

controllers. Assuming this and the satisfaction of Eqs. (51,52), becomes 

V a 7 a X0 a X8° 6 + 2la X a 3 a X0 ' a 3 a 7 ' “xV^ ’ (59) 

Hence, for non-zero values of fl, small , and with the geometry given by 
Eqs. (55,56), the condition for non-negative values of can be obtained as 

4T^ + 9M 11 M 22 B 2 >. 0 , (60) 

The above relation is satisfied almost certainly for all real satellites, 

A similar treatment for the coefficient yields the inequality 

X 4 

I* - 2M 11 M 22 n 2 - — [(Xf^ 2 - (Xf 2 ) 2 ](M 23 -M 13 ) > 0 (6X) 

which is also satisfied. 

Carrying on with this procedure, it can be shown that the coefficients 
a 4* a 3 * a 2> an< * “l are well behaved and positive definite. Thus, the 

only critical coefficient is oig. This is approximately given by 

a 8 " ~ a 10 n m ~ ZI ^ M ll +M 22^ (62) 

x 

where and are the variances of the inertia noise along and 1^, 
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respectively. 

To give a clearer picture of the error growth phenomenon , we will analyze 
the response of a three-axes stabilized satellite. 


Let it be assumed that initially 

fl » 0 
a c * 0 


Xf 3 - e 6 - 0 


(63) 

(64) 

(65) 


and Eqs. (51,52) are satisfied. In this case, all coupling in Eqs. (32-40) 
are lost and the responses grow linearly with time, according to the relations 

w - « 


**100 “ 1**1 “ a 2 ^ 

M 010 * [ * F 2 ~ a 4 ]t 

,2.2 


200 


M. 


020 


[Af 1 - a 2 ] t 
[Xf 2 - a 4 ] 2 t 2 


( 66 ) 


M 002 * 0 

The growth rate of the responses is greatly changed if Eq. (65) is not used, a 
though Eqs. (63,64) and Eqs. (51,52) are used. In this case, the following 
four equations remain coupled: 

*100 “ a l *011 

M 011 “ * Xf 3 -a 6**010 + ^*2 _f V*001 


where 


M 010 “ ' a l M 10l 


(67) 


M 101 ” ^ X *3 -a 6^ M 100 + ^l -a 2* M 001 


M 001 ’ (Xf 3- a 6> 


The eigenvalues of Eq. (67) satisfy the following algebraic equation; 
[p* * a 2 (Af 3 -a 6 ) 2 ] - 0 


( 68 ) 
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Equation (68) states that, apart from the linearly growing components, there 
will be exponential and sinusoidal components in the satellite response, when 
(XF 3 -a 6 ) is large. 

The above mentioned cases, identified by Eqs. (66) and (67), are extremes. 
A real situation can be portrayed better by assuming (XT^-a^) is non-zero but 

Very small, leading to a slight coupling in Eqs. (32-40). This causes a small 
non-zero value of ft to be developed, although Eqs. (51,52) are satisfied. With 
this compromise, the eigenvalues of Eqs. (32-40) satisfy the following charac- 
teristic equation: 

i p 2 (p 2 +fi 2 ) 2 [p 3 -a 10 n 2 p 2 + (4n 2 -a 7 a 18 n 4 )p + (a 7 a 10 a 18 n 2 -2a 7 -2a 18 )n 4 ] 

- a 5 Kp - 0 (69) 

where it la assumed that 1^ “ I 2 “ -j Ij and 

K - (X 2 -X 2 )p 4 - 8X 1 X 2 «p 3 + p 2 n 2 (tt 2 -A 2 )(5+a 7 a 18 n 2 ) 

+ 2A 1 ?i 2 n(a 7 +a 18 )] + ApnV^d+a^gO 2 ) + 

f + 0 5 [a 7 a 18 fl(^-A 2 ) - 2A 1 A 2 (a 7 +a 18 )J (70) 


I In Eq. (70), X and X 0 are given by 

I 1 

\ \ - ^ - a 2 (71) 

} 

I * 2 ” ” a A • C72) 

Equation (69) can be viewed with a better perspective by considering a_, 
a^, and a^g to be small. This reduced Eq. (69) to the form J 

p 3 (p 2 +SJ 2 ) 2 (p 2 -a 10 n 2 p+4!J 2 ) - 0 (73) 


| It is now clear that a spinning satellite will begin to satisfy Eq. (73) 

j immediately In the presence of noise. A three-axes stabilized satellite, on 

» the other hand, will satisfy Kq* (73) only after ft period of linear error growth. 

If a 0 is equal to zero, Eq. (73) predicts a dominant cyclic response with the 
well 1U known frequencies of n and 20 . The solutions of Eqs. (32-40), corres- 
ponding to the characteristic Eq, (73) , are easily obtained as follows: 


M 001 

**002 


- [J 





17 


where 


M 100 - - jj- + (A 1 +A 3 t)cosnt + (A 2 -A 4 t)sinftt 
Mq^q ■ "jj" 4* (A^+Ajt^slnflt “ (A 2 “A^t) cosflt 


M 101 * AjSlnOt + A^cosGt 


M 011 * A 4 sinnt “ A 3 cb»nt 

M 20 o * c + exptj a 10 fi 2 t](D^c682nt + D 2 sin2Gt} + G^(t) 

Mq 2 o “ C - exp[|a 10 n 2 t]{D 1 cos2«t + D 2 sin2f2t} 4* G 2 <t) 


(74) 


110 


A-X exp[i a 10 fl 2 t] 

+ ■ ~ {( 16 D l^a 10 «D 2 )sln2Gt 

G (l6+a 1() G ) 


- (16D 2 +4a 1() GD 1 )co82fit} + Gj.Ct) 

, A 2 , A^, A^, C, and D 2 are arbitrary constants, and 


G^t) 


a ift X l X 2 ^1 

— 2U " + V*!! 8 * 11 ** 4* r 12 cosGt] 

+ A 2 [r 2 ^sinGt + r^cosGt] + A^tr^tsinflt + r^tcosGt 
+ r 33 sinfit + r^cosQt] + A^fr^tsinGt + r^tcosflt 
+ r^sinGt + r^cosfit] (74a) 


G 2 (t) 


a l0 X l X 2 X 2 

- — - 4- A 1 ts 11 cosGt 4- s 12 sinGt] 


2G 


+ A 2 [s 21 cosGt 4- 8 22 sinftt] 4- A^fs^^tcosGt 4- s^tsinGt 
4* s^cosftt 4* s^sinGt] 4- A^[e^tcosGt 4* s^tsinftt 
+ s^cosGt 4* s^sinfit] 

G 3 (t) - A^t^cosftt 4- t 12 sinGt] 4- A 2 [t 21 cosGt 4- t^sinGt] 

4- A 3 [t 31 tcosGt 4- t 32 tsin«t 4- t^cosGt 4* t^sinfit] 

4- A^[t^^tcosGt + t^tsinfit + t^cosftt 4- t^sinfttj 


(74b) 


(74c) 
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44 * 

[81X 2 « J + 

72X^7^ - 21X 2 fi7T^ - 

2 X^1 


’41 " 

“ 2t 41 

5 r 42 ‘ 

2X i 

' 2t 42 + [i 5 

r 43 " ‘ 

' 2t 43 " 

2t 42 

a 


2t 44 ' 

2t 41 





2X 2 

‘44 * 

n 

; s 4i * 

- 2t 42 1 

; s 42 ‘ 

Zt 41 + 

n 

43 * 

-2t 44 


+ i. 
« 2 ■ 

1 s 44 “ 

2t 43 + 

2t 42 

n 



(74h) 


The nature of the functions G^(t), G 2 (t), and G^(t) can be given a simpler form 
if a^g is neglected in Eqs. (74a-74h). In this case, the functions are given by 

* 

G l (t) “ " "I + 'T [X 1 A 3 - X 2 R(A l +A 3 t)]cosnt 

- \ [A 1 A 4 + x 2 n(A 2' A 4 t)lslnn,: 

X 2 2 

c 2 ( t) - - -f + -7 [* 2 a 4 " x 1 n(A 2 -A 4 t)]co 8 nt 

2 (74i) 

+ ~2 lA 2 A 3 + (Aj+Ajt) ]sin£Jt 

G 3 ( t) - [X 1 (A 1 +A 3 t) + X 2 (A 2 -A 4 t) + | (X 2 A 3 +X 1 A 4 )]co 8 nt 

- ^ [X 2 (A 1 +A 3 t) - X 1 (A 2 -A 4 t) - i(X 1 A 3 -X 2 A 4 )] 8 inflt. 

The constants A^, A^, A^» A^, C, and D 2 are calculated from the appropriate 

initial conditions. Equations ( 66 ) and (74) provide a basis for comparison of 
the error growth rate of spinning and three axis stabilized satellites. If a^g, 
given by Eq. (62), is large and A^ or X 2 are small, then a three-axes stabilized 

design is warranted. The reverse is also the case. Interestingly enough, all 
these predictions have been borne out by analog simulations. 


THE PERTURBATION SCHEME 

A perturbation solution of Eq. (1) will now be obtained with the assumption 
that A 1 ^ and that the noises associated with the moments of inertia of the 
satellite are small. The inertia noises are defined as 

I ± - T a + e ± + n t (t) ; i « 1,2,3 (75) 

where and are the noises in the sample and time spaces, respectively. 

The angular velocity responses, are assumed to be functions of the seven 
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small parameters A*, and of the form: 

“i “ a t + X '“i0 + E l"il + e 2“i2 + E 3 w i3 + n l“i4 + n 2“i5 + n 3 u i6 
+ (X’) 2 u i7 + X' Cl o> 18 + X*e 2 u 19 + X’ e 3°> il0 + 

+ V Vll2 + X'n 3 u il3 + (^) 2 u 114 + + Wilt 

+ c 1 n 1 n> il7 + Wil g + Wll9 + (e 2 ) 2 6. 120 + e 2 e 3 u i 21 

+ e 2 n 1 ii) i22 + e 2 1 2 ^ i23 + E 2 n 3“i24 + (e 3 } u i25 

+ E 3 T 'l“i26 + Wi27 + E 3 n 3 U i28 + (ll i )2u i29 + W±30 

+ , 'i T >3“i3i + (l 2 ) “ i32 + T '2 I, 3 u i33 + ^ n 3^ u 134 * 

In Eq. (76), the cubic and higher powers of the small parameters are neglected. 
The quantities ft* are the nominal values of the angular velocities ok. It is 
assumed that 

fij * n* « 0 


ft* * ft « a constant 

03 1 (0) - U) 2 (0) » 0 =* [o> 3 (0) - ft] 


(77) 


Equations (75), (76), and (77) are substituted into Eqs. (1) and separate 
equations are then formed corresponding to each of the various combinations 
of the small parameters. This classical principle of separation of parameters 
results in only a few of the multitude of terms on the right hand side of Eq. 
(76) being non-zero. Thus, a more compact expansion for the angular velocities 
is obtained as 

“x * X ’“l0 + (X '> 2 “17 + X ' E 1 U 18 + X ' e 2 U) 19 + X ’ E 3“ll0 + X 'Vlll 


+ x , n 2 u 112 + x'n 3 u 113 


2 

XVw 2 o + (A 1 ) w 2? + X , e 1 (n 28 + A 1 e 2 tJ 29 + X,Er 3 w 


210 


+ X,r) l W 211 + X **Z a HZ 4 X * ^3^213 
2 

w 3 * ^ ^ * W 3Q + w 37 + X * e 3 a> 3l0 + X * T1 3 U) 313 


Let L* be a matrix differential operator defined by 
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1 




— _ 


1 

T x 0 0 

1 

“1 


0 d 3 -i 2 )n 0 

1 

“1 

I L*(m^) * 

0 i 2 0 



+ 

dj-ipn 0 0 


“2 

1 

001. 
k ■* 

. 

1 

“3 

1 

O 

O 

O 

1 

1 

“3 


( 79 ) 


Then the perturbation equations for the components of jw given in Eq. (78) take 
the form 

L*(A'u i0 ) - [X 1 ^, X'f 2 , X'f 3 ] T (79) 

L*[(X') 2 ai i7 ] = [ (I 2 -X 3 ) (X ' ^ 20 ) (A ' w 3Q ) , (13-ip (X'iu 10 ) (X ' t^) , 

(Y 1 -I 2 )(X*(o l0 )(X'<« 2O )] T (80) 

L * (A ' e l“i8 ) " ^ t{X , E 1 0(T 3 -I 2 )u 20 - A'e^} , -X'e^SJai^, 0] T (81) 


2 u i9 ) 

= i (X' e.I.na: 
I 2 ^ i 

' 2 q> “A £^1 ( I 3“l-^)^ w io ^ ^ 2 } * 

(82) 

3“il0 ) 

= [-ne 3 x' W20 

, «te 3 x»*. 10 

’ " e 3 X ’“30 ]T 

(83) 


= - [X’n^io 

* An^io * 

oi 1 

(84) 

l 2“il2 ) 

* [A , n 2 fi w 20 , 

“^ ,n 2 to 20 f 

0] T 

(85) 

l 3“il3 ) 

= (—x'n^owgQ 

> A , n 3 «W 10 

, -A ! n 3 w 30 ] 

(86) 


Equations (79-86) are easily solved. In particular, assuming 

X - Y = i I 
1 x 2 2 3 


A'f 2 * X’f 3 » 0 


(87) 

( 88 ) 


the solutions to Eq. (79) and (80) are obtained as 
1 

A f u) in = — / cosfi(t-t)A , f (x)dx 
iU \ 0 - 1 


A*0J 9n = - — / sinfi(t-x)A , f 1 (x)dx 
X 1 0 1 


(89) 


30 


“l7 “ “27 


“37 = ° 


and hence 


“310 = “313 = 0 


(90) 
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The perturbation solutions obtained so far from Eqs. (89,90) agree closely 
with the Fokker-Planck solutions given by Eq. (74). But the drawbacks of the 
perturbation scheme become apparent when Eqs. (81-83) are solved. Equations 
(81-83) predict a secular growth of the angular velocities even for the time- 
independent sample space inertia noises, e,. This is obviously not true from 
a physical standpoint. Thus, all perturbation equations involving e., but not 
n ± , must be discarded and the parameters must be absorbed in 1^. equations 

(81-83), then, are discarded and are set equal to zero, so that Eq. (78) re- 
duces to 


'“lO + A 'Vin + v 


" 2*112 + X ' n 3“ll3 


X ’“ 2 0 + X 1 


V211 + X 


" 2*212 * X n 3 U 213 


“3 = 55 + X ' n 3“313 


(90a) 


Equation (90a) predicts that, if A’f^ and are independent, then the mean val- 
ues of the amplitudes of u) 1 and do not grow with time. It also states that 

the variances of the amplitudes are stable and oscillatory and that the ampli- 
tudes of oscillation of the variances are constants for all time. In other 
words, no growth rate of the variances of is predicted by Eq, (90a), Con- 
trary to this prediction, it will be seen in analog simulations that the ampli- 
tudes do grow with time, even if and r } are independent. 


THE MOST-LIKELIHOOD APPROACH 


The method of most- likelihood estimates will now be applied to the system 
described by Eq. (1). As mentioned earlier, this method is based on maximizing 
the joint probability density of the random variables under the constraint 
that Eq. (1) holds. It can be shown that this method, when applied on even a 
linear equation, finally requires the solving of a nonlinear equation. For 
this reason, the nonlinear Eq.. (1) needs to be linearized Initially to make ana- 
lytic manipulations possible. 

The well-known linearized form of Eq. (1) is given by 


where 


Vi 

z 2*2 


h*3 



= \'f 1 - (I 3 - I 2 )Sto 2 

f 2 = x'f 2 - <i r i 3 )«ta 1 


( 91 ) 


( 92 ) 


Let v be the vector defined by 
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i 

j 


V - [v 1 ,ii 2 ,y 3 ,6 1 ,« 2 ,s 3 ] T 


(93) 


Let the matrix elements be defined by 





E iv ± Vj} - 1^6(0 


(94) 

. 

Let the functional J be defined by 





J “ i.j./o {V i lN ’ 1] iJ V J + 2P k [I kV F k 1)dt 

(95) 


. 

| where 3j are arbitrary time-dependent Lagrangian multipliers. It can be shown 

j UO] that the most likelihood estimates of u), are obtained by minimizing the 

' functional J in the interval [0,T] with respect to the variables and 


The variational equations for minimizing J 
following two equations: 

are given by Eq. 

(91) and the 



j tN "\j V j + 4k^ (I kVV 

= 0 

(96) 


■ 

[ • 

A 3F 1 

5‘W + j ° 


(97) 


j The terminal point condition on 0 is given by 




1 

Pjd) = o. 


(98) 


| Assuming that 





S 13 “ N 23 ” N 31 " N 32 * 0 


(99) 


and 

% ,= N j4 = 0 

N 5j ■= N j5 = 0 if j ji 5 

% “ N j6 “ 0 lf J * 6 


(100) 


Equation 

(96) can be opened up to read 





W 1 “ N 11 [f Vl + “*1^ + N ’l2 IB 2“2 " 





u 2 = N 12 [6 -jU^ + nw 1 $ 2 ] + N 22 [3 2 w 2 - 

P 3 * N 33 [^ 3 w 3 + 0(6^2 - ] 

S 1 - - N 44 e i 

06^2 ] 

(101) 

1; ' 


2 4 


S 2 “ - N 55 0 2 

{ 3“- N 66S 

Using Eqs. (92) and (101), Eqs. (91) and (97) can now be reduced to the follow- 
ing forms: 

[Ti + N 11 (e 1 m 1 +nm 1 B 2 ) + N 12 (e 2 iii 2 -!2B 1 o) 2 )]oj 1 + r[T 3 ~T 2 

+ N 33 C B 3 «>3 + (U&jWj-iDj^) ) - S 12 ( p 1 “l +nu 1 p 2 ) 

- N 22 (B 2 “ 2 -R6 1 u 2 ) ]u) 2 - X' f x + N 4 a B 1 = 0 (102) 

[1 2 + n 12 (B 1 “i + n “x e 2 > + N 22 (e 2“2" n6 l“2 ) “2 + nt *l“*3 

- VW^lV + N 12< 6 2“2“ nB l“2 ) - N 33 {B 3“3 

+ n(B 1 u 2 -0) 1 B 2 )}J(u 2 - Pf 2 + N 55 B 2 = 0 (103) 

[T 3 + n 33 {b 3 iL 3 + n(B 1 K 2 -u 1 B 2 )}]a! 3 - X'Tj + n 66 b 3 ■= o (104) 

[II + n 11 (B 1 “i+«“ 1 B 2 ) + N 12 (B 2 i 2 -!2B 1 m 2 )]B 1 - 

+ N 11 (8 1 m 1 +nm 1 B 2 ) + ^(B^-RB^) - 

+ f2(3 1 w 2 -uj 1 3 2 )}]3 2 « 0 (105) 

[ 1 2 ■+ N 12 (e A +0M l $ 2 ) + N 22 (g 2V fiP l a> 2 )] ^2 " nt V*2 

+ V e 3“3 + “ (P lV W l P 2 )} - N 12 (P A + ^1 P 2> 

" N 22 (B 2“2'" ne i ai 2 )l3 l “ ° (106) 

[X 3 +-»3 3 {&3« 3 . + fi(e i aJ 2 -u> 1 0 2 ) >]3 3 - 0 (107) 


Equations (102-107), together with the Initial conditions on oo^ and the end 
conditions on 8^ given by Eq. (98), form the final two-point boundary value 
problem Governing the stochastic motion of the satellite. To solve this 
problem, a perturbation sequence for 8^ and has to be adopted. 

Let it be assumed that e is a small parameter and the numbers N,, are of 
the order of e or less. Let N^* i = 1-7, be defined as 3 

Nn * £N 1 

n 12 - eN 2 ■ 
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N 22 “ eN 3 
N 33 “ eN 4 
N 44 ‘ eN 5 
N 55 “ eN 6 
N 66 “ eN 7 


where 


Max |N 


ij 1 


Let the variables and 3^ be assumed in the form 


“i ■ “ i0 + E “il + E “12 + ••• 


B i“ B iO +E6 il + E6 i2 + •" 


such that 


(108) 


(108a) 


(109) 


B ±j (T) * 0 (110) 

Substituting Eqs. (108) nnd (109) in Eqs. (102-107) and separating the co- 
can be seen that the zeroth order response 

( 111 ) 


- A’f 2 (U2) 


After some involved algebra and the use of Eq. (110), it can be seen that the 
predicted response from the higher order perturbation equations has essential- 
ly the same characteristics as that obtained by the straight forward pertur- 
bation scheme explained in the preceding section. Thus, the method of the 
most likelihood estimates suffers from the same drawbacks as those of the per- 
turbation method. 


THE METHOD OF STOCHASTIC EIGENVALUES 

According to this method, the deterministic solutions of Eq. (1) are to 
be obtained first. Stochasticity is then imposed on these solutions to esti- 


efficients of e , e , e , etc., it 
is given by 


B i0 “ 0 

I l 6? 10 + ^ I 3’ I 2^ fiw 20 


I 2 U 20 ” ^ I 3“ I l^ a) lO 


V30 " X?f 3 
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mate the behavior of the system which was random from the beginning. Although 
this method is not exact, it is much simpler than the methods previously dis- 
cussed. 


?or example, the approximate deterministic response of a three-axes stabi- 
lized satellite is given by 


t A'f 

0 ). « / dt , i* 1,2,3 
1 0 L ± 


(113) 


Hence, assuming to be a constant, the mean values and the variances of 

are given by [16 ] 


where 


Vf 

t 


E{oj^} * (a ± ) 2 t 2 


(114) 

(115) 


2 1 r 

°l “=2 l ~ 


(X’f^ 2 ^ + X^M 




4 + M n 


2 1 


(X'f^ 2 ^ + I ^ 1 55 , 
*■ —2 
4 + «22 


(116) 


2 1 
a 3 = —2 

4 


a , f 3 ) 2 M33 + I^M 6 

*■ O J 


I 3 + M 33 


In deriving Eq. (116), it was assumed that and 6 ^ are Gaussian random variables 

For the case of a spinning satellite with 1^ = ^ c y *^3 * anc * 
constant values of T'T 1 and T'f^, the deterministic amplitudes and frequency of 
oscillation of and oy are given by 


Freq. [u^] 
Amp. [u^] 

Amp . [a> 2 I 


‘ Freq . [a) 2 3 = 

_ *' f 2 
I 2 m 3 r 2 n 
X'fj x'f x 
*X*3 T l n 


(117) 
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when m^(0) * u> 2 (0) - 0. 

Hence* the growth rates of the amplitudes and frequency are described by 
the variances, which are 

EtlFreq.^]] 2 } » E{[Freq, [u 2 ]] 2 > - (o 3 ) 2 t 2 (118) 


, (X'f ) 2 o 2 t 2 + a\. 

E{ [Amp . [to.]) 2 } - — 1 

i 2 n or + <j 2 t 2 

2 1 < X ' T l) 2 °3 t2 + a \t, 

E{[Amp. [uj]] 2 } - rfr 


i 2 n 2 


n 2 + o 2 t 2 


(119) 


From Eqs. (114) and (115), it is seen that the approximate predictions 
for the responses of three- axes stabilized satellites are quite satisfactory. 
Equation (118) approximately predicts the frequency growth phenomenon. Equa- 
tion (119) predicts that, when t is small, such that o^t is small compared to 

ft, the variances are of the form 

(X’f 2 ) 2 o 2 t 2 

EUAmp.[<o jr> - + M ] (120) 

1 i 2 n 2 fi 2 33 


But for large values of t, the variances will reach a constant value. This is 
given by 


EUtap.i^n 2 ) = (^|) 2 


( 121 ) 


The prediction of an initially growing variance finally levelling off to 
a constant value is satisfactory and is corroborated by analog simulations. 
The only problem with Eqs. (120) and (121) is that these equations predict a 
lower growth rate and a lower value of the asymptotic variance as ft becomes 
large. In this respect, Eqs. (120) and (121) differ from the Fokker -Planck 
formulation and the analog simulations which give higher growth rates and 
higher values of the asymptotic variance for larger values of 0. 


ANALOG SIMULATION 

The results of simulation of the satellite response, as given by Eq, (1), 
can now be presented. The simplified system block diagram is shown in Figure 
1. This system was programmed on an AD-256 (Analytical Dynamics-256) analog 
computer. The white noise Inputs p. and , i « 1-3, were obtained from a 
coupled SDS-930 (Scientific Data Systems-930) real time digital computer. A 
high frequency RO (Repetitive Operation) clock circuit from the AD-256 was 
used to trigger a pseudo-random number generating program in the SDS-930, Sam- 
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pies of twenty such pseudo-random numbers were used to form a Gaussian white f 

noise sequence with a zero mean value and suitable peak values. Six such in- 1 

dependent noise sequences were continuously generated in the 3DS-930 and fed j 

to the AD-256 through six DAC (Digital to Analog Converter) lines. One test J 

line was also used to interrupt the SDS-930 and change the peak values of the I 

noise sequences. A sample of the noise sequences p , i * 1-3, is shown in | 

Figure 2 at a high brush recorder speed. At any instant of time, the frequencies f 

of generation and the peak values and, hence, the bandwidth of all p^ and 6^ | 

i * 1-3, were maintained equal. Thus, 6 , i - 1-3, are similar in nature to | 

that shown in Figure 2, although all six noise sequences were independent of *} 

each other. | 

jj 

Let r be the ratio defined by r 

r « [Peak value of p and 6 , i * 1-3] /I^ (122) 

114 1 

where I 3 is the nominal moment of inertia about the spin-axis. Brush records fj 

of the simulated angular velocities uv^, and for different values of fj 

r and ft, are shown in Figures 3-15. The values of r and ft, corresponding to | 

each of these figures, are tabulated in Table 1. In all cases the initial | 

values of and were taken to be zero. | 


Table 1: Index to the attached figures showing samples of the sto- 

chastic satellite responses. 



EVALUATIONS AND COMPARISON 


The results of the analog simulation will now be evaluated and compared 
with the predictions of the analytical methods discussed earlier. 


The first important result of the simulation study is that, in every case, 
the responses grow with time. The growth phenomenon is predicted by all four 
of the analytic methods only for the case of a three-axes stabilized satellite. 
This was true because, if ft » 0, Eq. (1) leads to a perturbed equation given by 















29 


Vi 


; ± = 1,2,3 


(123) 


Responses given by the solutions of Eq. (123) are the integrals of X , f. and, 
hence, must grow linearly with time, even if X'T^ are equal to zero. nut, 
in the case of spinning satellites, only the Fokker-Planck formulation pre- 
dicts an initial exponential growth. The perturbation method and the most- 
likelihood approach predict a constant variance. The stochastic eigenvalue 
method also predicts a linear growth rate which, however, is inversely propor- 
tional to ft2. Looking at Figures 4, 8, and 12, or at Figures 5,9, and 13, 
or at Figures 7, 11 and 15, it is seen that the variances increase with ft. 
Thus, at this point, the Fokker-Planck formulation is apparently the best of 
the theories under consideration. 


A second interesting result, discernible from Figures 3, 7 and 11, is 
that, with time, the response amplitudes reach a stable value. Such stable 
values are predicted directly by the stochastic eigenvalue method. The per- 
turbation method and the most-likelihood approach also yield the same result 
if it is assumed that these methods are valid only for the asymptotic case. 

It is to be noted that the Fokker-Planck formulation can also be made to 
yield this result, although not as directly as the other methods. To do this, 
let the solutions of ^CK) an< * **020 aS 8 * ven ^ '(^4) be considered: 


M 200 = C + exp[-| a 10 n Z t]{D 1 cos2nt + D 2 sin2f!t} + G^t) 
Mq 2 o “ C - exp[-| a 1() n 2 t]{D 1 cos2nt + D 2 sin2fit} + G 2 (t) 


(124) 


The exponential terms in and M Q2 q appear with opposite signs. 

According to Eq. (124), one of the variances must grow and the other de- 
cay with time. Thus after a certain time, one of these variances will tend 
to be negative. But variances are by definition non-negative quantities. 
Hence, and D 2 are to be taken as non-zero until one of the variances first 

becomes zero. and should then be set equal to zero in order not to have 

negative values of ^ 2 Q0 anci M 020‘ This P roce<3ure yields the prediction that 

the response amplitudes become stable after a certain time, which is in agree- 
ment with the simulation results. 


The last obvious result obtained from the simulation is that, for a given 
value of ft, the variances and the growth rates Increase with r. This is ex- 
pected , both intuitively and rationally, and all four theories predict it. 

A comparison can now be made of the theoretical methods, based on purely 
analytical grounds. The strength of the Fokker-Planck method lies in the fact 
that it does not require either uncoupling or linearization of coupled non- 
linear systems such as that of Eq. (1). The statistical moments of all orders 
are obtained directly as the solution of a coupled linear set of equations. 

Hence, digital computer methods can be used easily to solve such equations. 

The other three methods are based on initial linearization and possible un- 
coupling. This linearization results in a loss of useful statistical information 
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There are, however, some disadvantages of the Fokker -Planck method. The 
primary disadvantage Is that all statistical moments are coupled. Hence, when 
the number of dependent variables is large, the resulting set of equations is 
more so, even if the third and higher order moments are neglected. This method 
then requires some foreknowledge of the higher order moments and the statistical 
forms of the input random functions. 

In view of the above discussion, the following conclusions can be made: 

i) The Fokker-Planck formulation yields the most complete information 
on the responses of a satellite with random disturbing torques and stochastic 
moments of inertia . 

ii) For a satellite with very small inertia noises, the spinning configu- 
ration is better than a three-axes stabilized configuration. The reverse is 
also the case. 

iii) In all cases, the responses have an initial fast rate of growth. 

But after some time, this growth rate falls off, leading to a constant variance 
level depending on the variances of the input disturbing torque and on the mean 
moments of inertia of the satellite. 
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STABILITY AND OONTBOL OF FLEXIBLE SATELLITES 


PART I - STABILITY 


T. C. Huang and Aniruddha Das 


ABSTRACT 

This investigation has two distinct parts . In this first part the environ- 
mental and control torques experienced by a satellite are assumed to be random 
so as to account for the Inherent errors in the control systems and the lack 
of exact models of the environmental torques . It has been shown that under this 
assumption the required stability criteria of a satellite is quite different 
from that obtained by a deterministic approach. It has also been shown that 
a flexible three-axes stabilized satellite can be made almost certainly asympto- 
tically stable, while the same is not true for a flexible spinning satellite. 


NOMENCLATURE 


A* « Composite body of a flexible satellite. 

[A. ] , i - 1-5 * Matrices associated with the equations of motion of the flex- 
1 ible elements; Eqs. (3), (49), (53) - (57). 

[AjJ, i * 1-5 - Matrices similar to [A^]; Eq. (44). 

a - Radius of the cylindrical rigid core of the assumed satellite 

configuration; Fig. 2. 

a* * Normalizing factor of the joint probability density; Eq, (17). 

B* * Additional composite body for a flexible dual-spin satellite. 


i\b [ b 2 3 

V b 3 
[Cl 

c ± , i - 1-10 


'**■ Matrices associated with combined equations of notion of the 
satellite; Eqs. (5) - (7). 

- Elements of r^ , Eq- (38) . 

* Stochastic system, matrix; Eqs. (21), (27). 

“ Coefficients of the characteristic Eq. (71); Eqs. (74) - (77), 
(80), (85), (86). 


e 1# 1 - 1-4 
K t) 



I*(t) 


Elements of r, , , i * 1-4 ; Eq . (38) . . . 

— ai 

Deterministic forcing function; Eqs. (5), (10). 

Elements of r^, r^; Eq. ( A) . 

Deterministic environmental torque vector on the satellite; 
Eqs. (4), (45), (50). 
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tG] 

g ij 

h [x(t) ] 
[I] 

[I] 


* Stochastic control matrix, Eqs. (21), (26). 

= Elements of coefficient matrix defined by Eq. (63). 

« Deterministic observed function of x(t) ; Eq . (13). 

■ Identity matrix. 

* Moment of inertia matrix of the nominal configuration of the 
satellite. 

A 

= Diagonal elements of [I]; Eq. (52). 

= The joints probability density of ( z- z) , (u- u) * (f^f) and 
[x (0) - x(0) 1 t Eq. (17). 


J* 

J** 

V i - 

[N ± ] , i “ 1-4 
to] 

‘V 

[P i ], i rn 1-5 

[Pj] ,1 = 1-5 
p i 

P*, 1 = 1-4 

[Q] 

Q ir 

a. 

a.' 

1 ■ x - 4 

[R] 


c Functional defined by Eq. (17a). 

= Functional defined by Eq. (18). 

= Lengths of flexible beams of the satellite. 

* Submatrices of [B^] Eqs. (28), (29). 

= Null matrix. 

■ Covariance matrix of [x(0) - x(0) ] ; Eq. (12). 

*= Matrices associated with the angular momentum equations of the 
flexible satellite; Eqs. (4), (50). 

* Matrices similar to [P^]; Eq. (45). 

* Eigenvalues of 

* Exponents of the assumed beam displacement function; Eq. (42). 
« Covariance matrix of [u(t) - _u(t) ] ; Eq. (15). 

* Elements of the characteristic matrix of Eqs. (49), (50); 

Eqs. (70), (71). 

* Generalized position vector of the flexible elements of the 
satellite; Eqs. (3), (4), (49), (50). 

* Vector, similar to Eqs. (44), (45). 

* Time dependent part of Eq. (42). 

* Covariance matrix of t_z(t) - ^(t) ] ; Eq* (14). 

* Displacement vector of the center of mass of the flexible sa- 
tellite from its nominal position; Eq. (63). 
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jr ,. , i “ 1-4 * Nominal position vectors of the spring-mass-damper systems; 

^ Eq. (38). 

, i * 1-4 “ Nominal position vectors of the beam-end masses; Eq . (37). 

[S K ], K = 1-3 * Coefficient matrix; Eq. (67). 

S = Covariance matrix of [f(t) - f(t)]; Eq. (16). 


T 

T 

T 

t 

_u 

u* 


X 

X 

*bi’ 1 ■ 1-4 

ydi* 1 = 1-4 

^di’ 1 “ 1-4 
1 “ x ' 4 

z 


JS 

a 

[a k ], k * 0-3 
6(t) 


1 

ii 



Generalized velocity of the flexible elements; Eq. (8), (11). 

Terminal point of controlling time interval. 

Terminal point of the time interval in which the maximum like- 
lihood estimates are required. 

Total kinetic energy of the flexible satellite. 

Time. 

Augmented control torque vector; Eqs. (5), (9). 

Control torque vector; Eqs. (4), (45), (50). 

Stochastic system state variable; Eqs. (21), (24). 
Deterministic system state variable; Eqs. (5), (8). 
Displacement vector of beams . 

Components of Eq. (41). 

Displacement vector of spring-mass- damper systems. 

Displacement vector of beam-end masses. 

Stochastic forcing function; Eqs. (21), (25). 

Observed values of the state variables; Eqs. (13), (19), (23). 
Characteristic values of Eqs. (49), (50); Eq. (70). 
Coefficients of structural equations; Eqs. (64) - (68). 

Dirac f s delta function. 

Relative angular displacement vector of A* with respect to B*. 
Lagrangian multiplier and state variable; Eqs. (18), (22). 
Lagrangian multiplier; Eq. (18). 

Coefficients of structural equations; Eqs. (64) - (66) . 
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T biJ 

• Beam displacement mode parameters; Eqs , (46) , (47) . 

T dl 

• Spring- mass -damper displacement mode parameters ; Eq 

a 

-Nominal angular velocity vector of the satellite. 

u 

- Perturbed angular velocity vector; Eq. (2). 

0)* 

-Angular velocity vector of the satellite. 

K 

- Angular velocity vector of A*. 

— B 

- Angular velocity vector of B*. 

OPERATORS 


C) 

» Time derivative; “■ . 

[ f 

■ Transpose. 

(') 

- Vector cross product operator; Eq. (69). 

o 

- Mean value. 

o 

« Vector. 

Det. [ ] 

- Determinant of the matrix. 

E[ ] 

« Statistical expectation. 

Tr. [ ] 

- Trace of the matrix. 

INTRODUCTION 



The primary requirement of an artificial satellite is that it should be 
capable of precise orientation in space. This capability is determined mainly 
by the stability and controllability of the satellite when viewed as a dynamic 
system. A large number of investigations have been made in the area of flexible 
satellite dynamics. But several interesting questions on the stability and 
controllability of flexible satellites have not been examined in sufficient de- 
tail. The present study looks at two of these questions: 

(a) What are the stability criteria of flexible satellites in the 
presence of errors in the controlling torques and largely unknown 
environmental torques? 

(b) For a given control system, and for a given number of torquing 
jets, is it possible to increase the controllability of a flexible 
satellite by monitoring the deflections of the flexible elements? 


In the first part of this study it will be shown that, in the presence of 
random errors in the external torques on a flexible satellite, the stability 
criteria are far more restrictive than those deduced from a deterministic 
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approach. The second part of this study will present reasons for an affirma- 
tive answer to question (b) • 

As mentioned earlier, deterministic criteria for the stability of flexible 
satellites have been studied extensively £1—4] . It must be noted that, to ac- 
count for errors in the external torques acting on the satellite, these torques 
and the dynamic state variables of the satellite model must be treated as 
stochastic variables. Several studies [5-7] on the state identification prob- 
lem have been done. These studies generally assumed Gaussian distributions 
and used Kalman filtering techniques. Using methods similar to that given in 
Ref. [8], equations of motion and the stochastic angular velocity response of 
flexible satellites have been computed in Refs. [9,10], But the problem of 
comparing the stability characteristics of various satellite configurations 
subjected to random excitations has not been investigated. 


DETERMINISTIC EQUATIONS OF MOTION 

Formal deterministic equations of motion of a flexible satellite can be 
established. The stochastic stability boundaries can be determined only when 
these equations are available. 

Let aj*(t) be the angular velocity vector of a flexible satellite. For a 
single body satellite, u*(t) Is a (3x1) vector. For a dual-spin satellite with 
two main composite bodies (A* and B*) , ^*(t) is usually taken as 

«*(t) - 0(t)] T (1) 

In the above equation, to* and jo* are the (3x1) angular velocity vectors of the 

composite bodies A* and B*; while _0 is the (3x1) relative angular velocity vec- 
tor of the body A* with respect to B*. Let be the constant vector of the 
nominal values of to*(t) , such that the perturbing angular velocity vector o)(t) 
is defined by 

m.(t) ® J£*(t) - (2) 

Let the motions of the flexible elements of the satellite be represented by 
the generalized (nxl) position vector jq(t) . With these definitions, the equa- 
tions of motion of the flexible elements can be expressed in the following fora: 

[A^iiCt) + [A 2 (u,u,n,t) ]cL(t) + [A 3 (u,i,£, t)]a(t) 

= [A 4 ](j (t) + [A 5 (u,n,t)]w(t) (3) 

Similarly, the equations for the conservation of angular momentum of the compos- 
ite bodies of the satellite can be shown to be of the form: 

[p^aCt) + [P 2 (w,a),n, t) ]a(t) + [ p 3 ti^t > 3a(t) 

« [P 4 ]m(t) + [P 5 (w,fi,t)]w(t) + u*(t) + f*(t) 


( 4 ) 
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where ja* (t) and J[*(t) are the controlling and environmental torque vectors, re- 
spectively. 

Detailed methods of developing Eqs , (3,4) are given in Refs, [1-4] and es- 
pecially in Reis. [11,12]. Eqs. (3,4) provide the coop let e set of equations 
of motion of the flexible satellite, Equation (3) contains in* scalar equa- 
tions, such that the matrices [A. ] , [k ^ ] and [A~] are square. Equation (A) con- 
tains either three or nine equations depending on whether the satellite is of 
a single body or a dual-spin type. 

Equations (3,4) can be conbined in the form 


iB^Jx + [B 2 ]x - u(t) + f{t) 
where, defining [1] to be the identity matrix. 


(5) 


and 



Ib 2 ] 



Aj- — A^ — A^ 


0 -I 


0 


x(t) - [u(t), s_(t), £(t)] T 


u(t) * [u*(t) , 0., 0.] T 
l(t) - [f*(t) , 0 , 0] T 


£<t) * £(t) 


( 6 ) 


(7) 


( 8 ) 

(9) 

( 10 ) 

(ID 


Equation (5) is the required differential equation describing the determinis- 
tic motions of a flexible satellite. 


STOCHASTIC EQUATIONS OF M5TI0N 

The stochastic equations of motion of the flexible satellite will now be 
obtained following the method shown in Refs. [8,9]. 

Let it be assumed that the initial values , x(0) , have a Gaussian distribu- 
tion with a known mean value, T(0) , and a known covariance matrix, [Pq] , given by 

[P 0 ] - E{[x<0) - 2(0)1 [X(0) - T(0) ] T } (12) 

Here the operator E denotes statistical expectation. Let jc(t) be monitored on 
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the Earth by measuring a variable z/t) where the mean value, z(t) , of z_(t) is 


related to js(t) by 

I<t) - h [x(t) ] (13) 

Let it also be assumed that the variables z^(t) t u(t) , and f^t) are Gaussian 
with known mean values and covariance matrices R(t) , Q(t), and S (t) , respec- 
tively. Hence, assuming zero lag, we get 

E{[z(t) - z(t)][z(x) - T(t)] T } - R(t)6(t-x) (14) 

E{[u(t) - u(t)][u(x) - t[(x)] T } - Q(t)S(t-x) (15) 

E{ (£(t) - T(t) ] [f (x) - T(x) ] T ) - S (t)fi (t-x) (16) 


where ^i(t) and f^ (t) are the mean values of u(t) and _f (t), respectively. 

Let the maximum-likelihood estimates of the response of the satellite be 
required in the time interval [0,T]. In view of the definitions given above, 
the joint probability density , J, of (z^ ~z) , (u-u) , (f-7) and lx(0) - x(0) ] is 
given by 

J * a*[exp(- —) ] (17) 

where J* is defined as 

J* - [x(0) - x(0) ] T [P 0 ] _1 [x(0) - x(0) ] 

+ /q [z(t) - I(t)] T [R(t)r 1 tz(t) - f(t)J 
+ [u(t) - u(t)] T [Q(t)] -1 [u(t) - u(t)j (17a) 

+ [f(t) - f(t)] T [S(t)r 1 [f(t) - I(t)] dt 
and * a** is the normalizing factor. 

The maximum- likelihood estimates can be obtained by maximizing the proba- 
bility density J. In other words, we minimize the functional J*, subject to 
the constraints that Eqs. (5) , (13) be satisfied. This is done by defining 
J** by the relation 

J** * J* + 2 /J “ hXx)J 

+ [ x + B^^{B 2 X-u-f_}]}dt (18) 

and minimizing J** by considering x(0) , z(t) , u(t) , f (t) , x(t) and the La- 
grangian vector multipliers ji(t) and _A(t) as the independent variables. 

It will now be assumed that 
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z(t) - h[x(t) 3 - x(t) (19) 

which means 

21 - [I] (20) 

9* 

With this assuaption, the variational, equations obtained by minimizing J** are 
expressed as 

X(t) - JC)X(t) + J6]u(t) + Z(t) , (21) 


X.(T) - 0 (22) 

and 

x(0) - x(0) + tP 0 ]V(0) (23) 

where 


X(t) 

Z(t) 


tx(t), X(t)) T 

IB" 1 f(t), - R -1 z (t) ] 



[C] 




B 


-1 

1 


[QfSHB' 1 ] 1 


(24) 

(25) 

(26) 
(27) 


Equations (21-23) are the required stochastic differential equations of motion 
of the flexible satellite. 


STABILITY CRITERIA 

The stochastic Eq. (21) has twice many scalar equations as the deter- 
ministic Eq. (5) . The deterministic equations are stable if the eigenvalues 

of have negative real parts. The stochastic equations are stable if 

all the eigenvalues of [C] have negative real parts. If there were no errors 
involves with u(t) and f(t), the matrices [Q] and [S] would be null matrices. 
Consequently , Eq. (21) would degenerate into Eq. (5). 

The hypothesis of this study is that [Q] and [S] are not null matrices, 
but have positive elements which are very small compared to those of (B,) or 
tBjIv Hence, half of the eigenvalues of tC] will be almost equal to the eigen- 
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values, p^, of [-B^ and the other half will he almost equal to -p^, That 

the eigenvalues of [C] lie symmetrically about the imaginary axis can be veri- 
fied by noting that 


TrlC] 


and that the eigenvalues of [B^Bg.]^ are a( 5 ua l an <* opposite to those of 

r -i 




In view of this, it is evident that Eq. (21) is always unstable. Even if 
the real parts of p. are zero, the instability will be caused by the multiple 
roots. Thus, according to the classical meaning of the term, no stability cri- 
terion exists for the stochastic Eq. (21). The physical reason behind this is 
that the probable errors in the dependent variables accumulate with time. This 
accumulation causes the maximum- likelihood estimates to be asymptotically di- 
vergent, even if the deterministic Eq. (5) is stable. The growth phenomenon, 
for a satellite in which the vector x.(t) is measured at discrete intervals of 
time, is illustrated in Figure 1. Let the mean values of ac(t) be considered 
to be given by the solutions of Eq. (5) . Let the variances of x(t) be compu- 
ted from the differences of the values of x:(t) computed from Eqs. (5) and (21). 
The error functions computed from these mean values and variances are shown at 
three instants of time in Figure 1. In Figure 1A, there is a data input and 
the computation cycle has been started. Hence the error distribution curve 
has a high peak. The variances here correspond only to the measurement errors 
of the variables _x(t) . In Figures IB and 1C, it is seen that the height of the 
error function becomes shorter and shorter, although the mean position given 
by Eq. (5) approaches the origin. In Figure 1C, the error function is very 
flat just before the new data input. It becomes sharp again just after the new 
data input when a new computation cycle is started. 

Since Eq. (21) is necessarily unstable, the stochastic stability cirteria 
for a flexible satellite must be formulated in a particular manner. The sto- 
chastic stability criteria of the response of a flexible satellite are those 
which make 

(a) the deterministic model given by Eq. (5) stable, and 

(b) the growth rate of the stochastic model given by Eq, (21) a minimum. 

In the absence of further information about the covariance matrices Q , R and S , 
these two requirements are met if the real parts of p^ are equal to zero. 

Thus, a flexible satellite will be called stochastically stable if all the eigen- 
values of [-B^ B^l ate purely imaginary. It is interesting to note at this 
point that a perfectly rigid satellite satisfies this requirement. 

Specific stability criteria can be obtained for a satellite when the ele- 
ments of [B^] and [B^] are known. For this, a particular satellite configura- 
tion has to be assumed. In the absence of such a specific configuration, sev- 
eral conditions sufficient to make the p^ purely imaginary can be established 
in terms of the matrices [A^] and [P^ ] , x i - 1-5, when the matrices are square. 
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SUFFICIENT CONDITIONS 

The sufficient condition for the p^ to be purely imaginary, the matrix 
must be an tisy metric. Let (A^ J and [P^], 1 * 1^5, be square 
matrices. Let [B~*] be given by 

tB x 3“ X 

Comparing Eqs. (6) and (28), the matrices N^, i - 1-4, are given by 


N 3 * 4 


0 
0 

0 I 


(28) 


iNj.1 ■ ip 4 - ?A\r l . 

[n 2 ] - (A 4 - A 1 P" 1 P 4 ]’ 1 

[n 3 ] - [V4V1!' 1 

in 4 i - [v 4 lp r A i 1_1 • 

Hence from Eqs. (7) and (28), (B~ 1 B 2 J ia given by 

- • *1 
[N 1 P 5 +N 2 A 5 ] -I n 1 P 2 +N 2 A 2 ] -t N 1 P 3 +N 2 A 3 1 

[n 3 p 5 +n 4 a 5 ] -[n 3 p 2 « 4 a 2 ] -{n 3 p 3 +n 4 a 3 1 

[0] [-1] [0] 

To have [b 7 A 1L] antisymnetric, the required conditions become 

h*S + N 2 A 5 “ 0 

n 3 p 2 + n 4 a 2 - o 

N 1 P 3 + N 2 A 3 ■ 0 
N 3 P 3 + N 4 A 3 ‘ ~ l 


[8- 1b 2 ] 


(29) 


(30) 


( 31 ) 


N 1 P 2 + N 2 A 2 * N 3 P 5 + *4 A 5 
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Eliminating i * 1-4 from Eqs. (29) and (31), the required sufficient con- 


ditions are finally obtained as 

[P 3 ] - tP 1 l (32) 

[A 3 ] - tA x 3 (33) 

[P 2 ] “ " (34) 

[Aj] - - [A 4 A^ 1 A 3 ] (35) 

[P 5 ] - (P 3 A" 1 A 5 ] (36) 


The stochastic stability criteria given by Eqs. (32-36) are much too re- 
strictive and it will be almost impossible to obtain a practical design of a 
satellite satisfying these constraints* For exanple, Eq. (33) requires that 
the natural frequencies of the flexible elements of the satellite should be 
equal to unity. This is not a feasible constraint. 

In spite of these drawbacks, Eqs. (32-36) do provide several guidelines 
for satellite design. It can be easily verified that Eqs. (34-36) are satisfied 
identically by a three- axes stabilized satellite in which all subbodies have 
undamped, purely elastic mountings. A spinning or a dual-spin satellite, even 
if it is free of damping, generally does not satisfy Eqs. (34-36). Equation 
(32) is satisfied by all types of satellites in which there is an axis of sym- 
metry, and in which the flexible elements are so constrained that the center 
of mass moves only along the axis of symmetry. Hence it can be claimed that, 
among satellite designs with comparable mass, stiffness, damping and covariance 
matrices, a symmetric, three-axes stabilized satellite is likely to have the 
lowest error growth rate. 


A SPECIFIC CONFIGURATION 

The constraints given by Eqs. (32-36) are too restrictive because, in their 
derivation, no attention has been paid to the Zero elements of the matrices in- 
volved. To utilize the location of the zero elements in the matrices [B^] and 
[B_] , a particular satellite configuration (shown in Figure 2) will now be con- 
sidered. The satellite consists of a rigid cylindrical body with four beams, 
four beam-tip masses, and four spring-mass-damper systems, placed symmetrically 
as required by Eq. (32). The beams are perpendicular to the axis of symmetry 
and are assumed to be axially rigid. The spring-mass-dat^er systems are assumed 
to be constrained to move only parallel to the axis of sym^try. These assump- 
tions lead to a large nunfcer of zeros in the matrices [B^] and [B^] , making the 
algebraic manipulations considerably simpler. 

The major drawback of any stability analysis with a particular satellite 
configuration is that conclusions drawn from it cannot be extended to other 
configurations. The method of modelling and analysis of the satellite configu- 
ration (shown in Figure 2) that has been used in this study partially overcomes 
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this disadvantage. In this method, the location of zeros in [B^j and [B^ 1 re- 
main unchanged when the nunbers of beams, tip-masses, or spring-mass- damper sys- 
tems are changed. 

THE DYNAMIC MODEL 

Let *a T be the radius of the main rigid body and 1 ■ 1-4, be the lengths 
of the beams.: Let and _r^, 1 * 1-4, be the nomLnalposition coordinates 

of the beam-tip masses and the spring-mass- danger systems, respectively. Ac- 
cording to the choice of coordinate axes shown in Figure 2, we have 

KaWj), 0, 0] T 
[0, -(a+O, of 

2 (37) 

[-(afi 3 ), 0, of 

[0, <a« A ), of 
Let it be defined that 

rji - tb x . 0, ei f 

^j2 m -fj* e 2^ 

(38) 

id3 “ t_b 3> °* e 3 ]T 
^4 “ I0 > V e 4 ]T 

Let x be the distance along the axes of the beams measured from the fixed ends. 
Let 2.^(0* and 9 * “ be tbe defiectlons of the beam-tip 

masses, the beams, and the sp zing- mass -damper systems, respectively. According 
to the previously assumed constraints, let it be defined that 

Z rl (t> = [0, y rl|2 (t), y rl> 3 (t)f 

*r2 (t) “ ty r2,l (t) * 0> y r2,3 (t)lT 

Z r 3 (t) - [0, y r 3 , 2 (t). y r 3 > 3 (t)f 

ix4 (t) ’ fy r4,l (t) * 0> y r4,3 (t)]T 


-rl “ 
-r2 * 
— r3 
-r4 “ 


( 39 ) 
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and 


- to. y bl ^ 2 (x,t), y blj 3 (x*t)j 
^(x.t) - ty b 2 > 1 (x,t), o, y b 2 ( 3 (x,t )] 7 
^(x.o - to, y b 3 > 2 (x,t), r^c *.*)! 1 
^ 4 (x,t) - ty b 4 > 1 (x.t), o, y Mi 3 (x,t)l 1 


i* ± (t) - to, o, y dl (t)] 


(40) 


(41) 


Equations of motion in the coordinates uj, y , . , y, . and y,. for i - 1-4, 

n»j Di,j oi 

j * 1,3 are obtained using the method shown in Ref. [1], The space dependence 
of these equations is eliminated by assuming 




[ q biJ ^) ] [ exp (P* x ) ~ Pjx-1] 


(42) 


and applying the Galerkin's method [1,11]. The space-dependent shape func- 
tions in Eq. (42) are assumed to be known and correspond to those of a canti- 
lever beam with a tip-mass. 


At this point, the boundary conditions 

y ri,j (t) ’ lexp(p *V ' P t *l " l fj 

are applied, and the equations of motion reduce to the form 

[AjJ<a/(t) + [Aj(w,m,n,t)]i , (t) + [A^(w,w,n,t) ]r’(t) 


[A^]w(t) + [A^(u,n)]u(t) 


and 


[Pj^ct) + [p 2 (]£>&> n,t)]^’ (t) + [P^(^,0),n,t)]^ , (t) 
* (p;]i(t) + [P|( 03 ,fi) Juj(t) + u*(t) + f*(t) 


(43) 


(44) 


(45) 


where jl* (t) consists of the non-zero elements of an< 3 y di > ^ ** The 

set of Eqs. (44) and (45) is of the order of 27. It is still quite difficult 
to extract any meaningful analytic stability criterion out of this set. 

It is now assumed that there exists certain unknown constants x,,. and 
i * 1-4, j ■ 1,3, such that 


x bl2 q bl,2 “ x b2i q b2,l ■ x b32 q b3,2 ” T b41 q b4,l 
T bl3 q bl,3 “ X b23 q b2,3 ” X b33 q b3,3 * X b43 q b4,3 


(46) 


( 47 ) 
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T dl y dl “ T d2 y d2 " T d3 y d3 " T d4 y d4 (48) 

The values of and can be obtained from the eigenvectors of Eqs. 

(44), (45). But it la not our intention at this point to look for eigenvalues 
and eigenvectors of Eqs. (44), (45). Substituting Eqs. (46), (47), and (48) 
into Eqs. (44) and (45), the equations of motion of the satellites are reduced 
to the form 


I A X^(t) + [A 2 li(t) + [A 3 J^(t) - l\]E + [A 5 1 w. 


(P^t) + [P 2 li(t) + [P 3 ] j l(t) « [P 4 ]«+ [P 5 )w +u*(t) + f*(t) (50) 


a(t) - [y dV % lf2 > %i >3 ] < 51) 

It should be noted that £(t) given by Eq, (51) is a (3x1) vector and all 
matrices [A^] and [Pjj, i * 1-5, are (3x3) matrices. The Eqs. (49), (50) now 

form only a ninth order set of ordinary differential equations. This great re- 
duction was made possible by the assumptions of Eqs. (46), (47), (48). It should 
also be noted that, irrespective of the number of beams or spring-mas s-daaper 
systems introduced at the initial stages of the dynamic modelling, Eqs. (49) 
and (50) can always be made a ninth order set by suitably augmenting the equa- 
tions in Eqs. (46), (47), and (48). 

Let it be assumed that the moment of inertia matrix, [I] , of the satellite 
is given by 


The linearized form of the matrices [A^3 and [P^, 1 * 1-5, can then be shown 
to be as follows: 


/ 1 * 
* w dl“ g 33 ll dl 


p bl2 -g ll lJ bl2“ g 21 ll bl2 ) 


^bl3 -g 32 w bl3* 

- g 32 w dl 
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U 2 ] 


[a 3 ] 


[a 4 1 


[A 5 ] 


[Pi] 


[p 2 ] 


-(8 ll W bl2 + 8 21 v bl2 ) 


0 


0 

2 

‘dl 


0 

0 

' J bl2 _8 ll* J bl2~ 8 

0 


0 

& 


0 

0 

0 

9 

y bl2 

d13 

v bl3 

0 

5 

'dl 

u dl 

0 

0 

0 

0 

‘bl3°3 

p bl3 n : 

3 0 

‘dlS 

"dA 

0 


p bl3 


( 54 ) 


( 55 ) 


( 56 ) 


(57) 


0 c* \ 

t 13 g 33“ S ll' ) 

0 

a 12 g 

( 0 4 . 

21 ^ a l3 g 32~ i3^ 

0 c* n 

*23 8 33 21 ; 

0 

a I2 g 

11 ^ a 23 g 32” S 23^ 

0 

, 0 ,0 
^ a 13 g ll +a 23 

g 2l“ S 32^ 

0 

.4 0 „ , 

, 2l“ a 23 g 33 ) 

0 

" a 12 8 ll 

^ S 23“ a 23 g 32^ 


l 13 g 33“ S ll 

0 

a 12 g 21 

/ 0 e 4 . 

(a 13 8 32" S 13 ) 

% 

0 

0 

0 



(58) 


( 59 ) 



K58ii!?»^ 
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( S 23 +a 23 8 33> 


“23 g 21 < S 32 +a 23 8 32> 


-(S 13 +o 13 g 33 ) 


^ S 3l"°13 8 32 ) n 3 


tP 4 J - [XJ 


[P,] - fltU - [in]. 


The undefined constants introduced in Eqs. (53) - (62) are defined by the fol- 
lowing relations: 


g ll 

0 

0 

q bl,2 ] 

g 2i 

0 

0 

q bl,3 | 

0 

g 32 

g 33 

y dl ) 


where _r is the displacement of the center of mass of the satellite from its 


nominal position, and 


P bl2 q bl2 + P bl2 q bl,2 “ y bl2 r c,l + v bl2 r c,2 + W bl2 r c,l 


+ W bl2 r c,2 + y bl2 r c,l + y bl2 r c,2 + w bl2“3 


^b!3 q bl|2 + W bl3 q bl,3 “ W bl3 r c,3 + y bl3“l + ,J bl3 u) 2 


+ y bl3“l“3 + y bl3“2“3 


V dl y dl + y dl y dl + y dl y dl “ P dl r c,3 + ^dl^l + W dl U 2 


+ y dl“l u '3 + ll dl UJ 2 <iJ 3 


1 O 1 1 

to ■ + r c,l [o ij ] ^ + r c ,2 [a ij + r c,3 [a ij ] ^ + y dl ls ij^ 


+ 1hl,2l S «lii + q hl.st S fi^ + 


T H bl,2 l ij ' M bl,3 l Ij T l “ij J c 

A 

where T is the kinetic energy functional [1] of the satellite. The operator 
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( ) on any (3x1) vector v Is defined by 

' 0 -v 3 v 2 

V - v, 0 -v. 


such that the cross-product between any two arbitrary vectors u and v is given 
by 

u x v * u v « -v u (69) 

Analytic search for the eigenvalues of Eqs. (49) and (50) is now quite 
easy, because these form only a ninth order set. As in the elements of the 
matrices [A^] and [P^], i * 1-4, these eigenvalues are functions of the unknown 

constants x^j and x^. The method of analysis to be adopted now is to obtain 
the stability criteria in terms of x^^ and x^. Then we must obtain the union 

of all criteria such that the resulting criteria become independent of x, , , 
and x di . b « 


EIGENVALUE EQUATIONS 

The characteristic equation in a for Eqs. (49) and (50) is given by 
[P 5 +oP 4 ] -[P 3 +oP 2 +a 2 P 1 ] 

Det. =0 (70 

[A 5 +aA^J -[A 3 +aA 2 +a A^] 


help of Eqs. (53) 

- (62) 

* it 

can be 

seen that Eq 



Q 12 

0 

Q i4 

<*15 

Q 16 


Q 21 

Q 22 

0 

^24 

Q 25 

^26 

Det. 

0 

0 

Q 33 

0 

Q 35 

0 

0 

0 

^43 

0 

^45 

0 


^51 

^52 

0 

<554 

0 

Q 56 


^61 

Q 62 

0 

<*64 

0 

Q 66 


It can be verified that the locations of the zeros of the matrix in Eq. 
(71) remain the same even if the number of beams or spring-mass-damper systems 
are increased. 
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Equation (71) can be factorized into 

'Qi 

Det. 

' Q 

Q, 

Thus, the characteristic equations become 


>11 

<>12 

Q 14 

<5X6 

>21 

q 22 

<>24 

<>26 

>51 

<>52 

«54 

<>56 

>61 

«62 

<>64 

<566 


IQ 43 Q 35 " Q 33 Q 45 ) “ 0 


and 


Det. 




(72) 


>11 

<>12 

<>14 

<>16 

>21 

q 22 

Q 24 

<>26 

>51 

<>52 

<>54 

<>56 

>61 

<562 

<>64 

<>66 


( 73 ) 


Equation (72) yields three roots of a and the other six roots are obtained from 
Eq. (73). One of the roots of a from Eq. (72) is identically equal to zero. 

The other two roots of Eq. (72) are given by the equation 


C i a + C 2 a + C 3 * 0 


(74) 


where 


C T * El z (y^^2”®ii ll bl2”®21^bl2^ + q8ii" a s*89i)] 


C 2 ” " I z (g ll U bl2 + g 21^bl2^ 
C 3 " I z (M bl2 -8 ll y bl2 _8 21 ll bl2 ) 


T»12 v 32 13rli "23 ■ 2l y 


(76) 

(77) 


Hence the requirement of purely imaginary roots leads to the conditions 

C 2 - 0 ; C 3 /C 1 > 0 if C 1 4 0 . (78) 

Expanding Eq. (73), the resulting equation in a is obtained as 

C 4 a 6 + C 5 a 5 + C 6 a 4 + C y a 3 + C g a 2 + C 9 a + C 1Q * 0 (79) 


To sitaplify the expressions of C^, 1 * 4-10, let it be assumed that 
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bl2 " ” T b32 ! T b21 ” " ‘b41 ’ ‘b!3 “ " % b33 * ^b23 " " ‘b43 1 


T v/.l * T ui o “ “* » T uo«a “ ” T v. 


T dl “ _ T d3 5 T d2 “ " T d4 


This assumed mode corresponds to that which, in terms of pointing accuracy, 
we are most interested. This mode leads to pure rotational motions of the rigid 
core about its center of mass, ffith this assumption, the coefficient C q is 
given by 


C 9 ’ <VV<VV 8 3 , ‘2l»il3 + 2b 1^13 n 3 S 23<W 

+ 2n 3( T d 2' ,T dl^ f 2 y bl3^ I x S 13 _S 21^ I z" I x^ ^ 

+ n 3 { S P bl3 V dl^ S 23 (I z _ V S 3lV ” I x S 31 (y bl3 U dl +f! 3 u bl3 W dl ) J 

- n 3 {y dl y bl3 S 13 (I z- I y ) + Ww'Was + .«4l<VY 1J 

+ S S 31 (I z-V ^ y dl y bl3 +n 3 V bl3 y dl^ ” S y dl W bl3^ I y S 32 +S 13 (I z _I y^ 

+ n 3 1 ^dl y bl3 +y bl3 y dr p bl3 y dl _1, bl3 ,J dl ) (S 23 S U _S 21 S 13 ) 

+ (y dl y bl3"' I bl3 W dl )(S 23 S 13 +S 32 S ll _S 13 S 23 +S 21 S 31 )1, (80) 

Expressions for the other coefficients in Eq. (79) are similarly obtained. 

For the roots of a in Eq. (79) to be purely imaginary, 

C 5 - C ? * C 9 = 0 (81) 

Examining Eq. (80) and similar expressions for and C_ it becomes evident 
that Eq. (81) can be satisfied for arbitrary values of 4^.. and if and 
only if ^ 

n 3 - 4 - 0 (82) 

Equation (82) is another proof of our previous claim that stochastic stability 
ii possible only for undamped three-axes stabilized satellites. 

THREE-AXES STABILIZED SATELLITES 

For a three-axes stabilized satellite, the constraints given by Eq. (78) 
are almost always satisfied. Also for this configuration, 


C 5 “ C 7 - C 8 ■ C 9 ‘ C lCf ° 
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such that Eq. (79) becomes 

or [Cg+C^a^] ■ 0 

Hence the required stability criteria are 


(84) 


0 1 c 4 


I X I y 1J dl v bl3 + ^VblsVtt + ‘‘JiVm + I y S ll p dl p bl3 


+ VdAl3 S n + < p dl p bl3 - 4 All' <4A3 ' s n s A < 85 > 


and 


0 - C 6 “ I y S 13 W dl w bl3 + I y S ll y dl y bl3 


+ I x S 23 W bl3 y dl 


+ 2b l y bl3 I x S 21 + I x I y (,1 dl 1J bl3 + * J dl P bl3 ) 


( 86 ) 


Constraints given by Eqs. (85) and (86) can be satisfied usually without great 
difficulty, irrespective of the values of and x^, This is due to the 

fact that and are mainly the mass and stiffness terms of the satellite 

model. Hence, it can be concluded that three-axes stabilized satellites are 
more likely to be stable under random environmental and control torques. 
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FIGURE I : INFORMATION STABILITY 



FIGURE 2- SATELLITE CONFIGURATION 


STABILITY AND OONTROL OF FLEXIBLE SATELLITES: 


PART II - CONTROL 


T. C. Huang and Aniruddha Das 


ABSTRACT 

This is the second part of an earlier investigation. In this section, it 
is demonstrated that, by monitoring the deformations of the flexible elements 
of a satellite, the effectiveness of the satellite control system can he in- 
creased considerably. A simple model of a flexible satellite had been analyzed 
in the first part of this work. The same model has been used here for digital 
computer simulations. 


NOMENCLATURE 
IA*K i - 1-5 


ftp, i 
ft('t) ] 

£ 

[X] 

CK] 

[K x l 


XA 


[ 0 ] 

[P p , i ~ 1-5 


- Matrices governing the equations of motion of flexible struc- 
tural elements of the satellite; Eq ♦ (1) 

= Matrices governing the satellite motion; Eqs . (3,7,8). 

= Upper (3x3) left comer submatrix of [$(t) ] [Bp 

« External forcing function; Eq : , (6) . 

- External torque vector on. the satellite; Eq . (2). 

= Identity matrix. 

= System fundamental matrix; Eq, (22). 

- Matrix defined by Eq. (28) * 

* Number of scalar elements In ^V. 

■ Null matrix. 

- Matrices governing the rotational motion of the satellite | 

Eq. (2) 


q\ {q[> 


- Generalized structural position coordinate vector, 
“ -Generalized position coordinate for the ith beam. 
* Terminal time for optimal control, 

.« Time, 
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’ T . 


t 
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U 

u, {u ± } 


u* 



X 


£ 

*ri 

A*, {A*} 
V*. (yJJ 

T 

m 

i\i> i - 

JW> (w i ) 

i f I, 
jd » IU 


* Torque magnitude parameter, defined by Eq. (25), 

- Generalized control vector; Eqs. (5), (22). 

* Control torque vector on the satellite; Eq. (2). 

- Various control torque functions; Eqs. (14) - (21). 

« State vector; Eqs. (3), (4), (22). 

* Uncontrolled response; Eq . (22). 

* Position vector of the ith spring-mass -damper system, 

* Position vector of the ith beam-end mass. 

= Control system parameter; Eqs. (24), (26). 

* Relative control torque magnitude vector; Eqs. (24), (27), 

* Dummy time variable. 

* Fundamental matrix of -[B|] E( l 8 • (9), (10), 

1-4 = Component matrices of [$>]; Eq. (11). 

* Angular velocity vector of the satellite; Eqs. (1), (2). 

* Various simulation responses of ij; Eqs, (14) - (21). 


INTRODUCTION 

In the first part [1] of this study, the question of stochastic stability 
of flexible satellites was discussed. Specific stability criteria were devel- 
oped for a simple flexible model of a satellite (shown in Figure 1) . In this 
part of the study, we determine whether it is possible to increase the pointing 
accuracy of a satellite by observing the deflections of the flexible elements. 
To do this, we use the same satellite configuration (Figure 1) and the theore- 
tical model developed in Ref. [1]. 

Likins and Fleischer [2] have shown that the flexible elements of space- 
craft can have a destabilizing Influence. They have shown a method of design- 
ing a proportional linear control system employing root-locus plots and eigen- 
value analyses. The control loop gains in [2] were based on a dynamic model, 
using hybrid coordinates, of a spacecraft containing long flexible beams. An 
essentially similar approach was employed by DiLorenzo and Santinelli [3], 

Here also a linear proportional control system was designed by considering the 
equations of motion of the spacecraft along with those of the flexible elements 
The spacecraft model in [3] consisted of a rigid body with two spring-mass 
systems . 

In this study , a time-optimal 1 bang bang* control policy has been assumed. 
The method of calculating the control torques is essentially the same as that 
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j; 


u 


r. 

jM 


given in Ref. [4]. Full details of the computation of control torques are pre- 
sented in Ref. [5]. Apart from the control policy, this analysis differs from 
Refs. [2,3] in another important aspect. In the analyses of Refs. [2,3], the 
deflections of the flexible elements are not observed. Hence, zero initial 
deflections and velocities of the flexible elements are inherently assumed. 

The present method can accommodate arbitrarily large initial conditions of the 
flexible elements of the satellite. 


THEORETICAL BASIS OF COMPARISON 

The theoretical analysis and comparison of the satellite responses is 
based on the dynamic model explained in Section 6 of Ref. [1], It was shown 
there [1] that, by using the Galerkin's method, the deflections of the flex- 
ible elements of the satellite are governed by purely time-dependent generalized 
position vectors, ^^(t), y ^ (t) and y^(t). It: was also shown that these vec- 
tors can be condensed subsequently , and reduced to a vector 3 / (t) by applying 
suitable boundary and continuity conditions. Usually the number of elements 
in J 3 L* is much smaller than that in the set 


Let a)(t) be the angular velocity vector of the satellite. Let Jj*(t) and 
_f*(t) be the control torque and environmental torque vector on the satellite. 
Given these definitions, it is well known [1,6,7] that the satellite response 
is governed by a pair of matrix equations of the form 

[A[]£'(t) + [A’(m,t)ji’(t) + [A^(m,t)]^(t) 

= [A’]m(t) + [A*(w)]w(t) (1) 

and 

[Pj&Ct) + [P^(w,t)]i'(t) + [P^(0),t)]^(t) 

“ [Pj]i(t) + [P * (w) ljw(t) 4- u*(t) + f*(t>. m 

Equation (1) governs the flexible motion of the beams, spring-mass-dampers, and 
beam-end masses of the satellite model. Equation (2) is based on the principle 
of conservation of angular momentum of the satellite. If (t) is a (nxl) vec- 
tor, then there are f n| scalar equations in Eq . (1). Equation (2) alwsys has 
three scalar equations. Equations (1) and (2) correspond to Kqs . (44) and (45) 
of Ref. [1]. 

Equations (1) and (2) are now combined together to form one first order 
equation given by 

[B^iCt) + [B|]x(t) - u(t) + f (t) O) 

where 

X= ,1 , ] T (4) 
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and 


u « [u*»0_,£]^ 





(5) 

f - (f*,0,0] T 





(6) 


[P 4 ] 

-[P{] 

to] 




- 

ia;i 

-IA^3 

10] 



(7) 


to] 

[0] 

II] 





"[p^] 

-IPJ] 

-IP!,] 



tB’] - 

[Aj] 

-[A'] 

-IAj] 


(8) 


[0] 

-[I] 

to] 




Let [ 4*(t) 3 be the fundamental matrix c 

>£ 

the homogeneous equation 


x - -[B| 


]* 




(9) 

that the solution of Eq. 

(3) is 

given by 



x(t) 


[♦(t)]x(0) + /{jl*<t“T)]lBjp~ 1 [u(t) + f(x)]dx . 


( 10 ) 


Let [ 4>( t> 1 be composed of [$^(t) ] , [^(t) ] , 311(5 t^Ct) ] such that 


[*] 


l\] [« 2 ] 

t» 3 ] 

(3x3) (3xn) 

(3xn) 

tV 


[2n x(2n+3) ) 



(ID 


when w(t) and jj/ ( t) are (3x1) and (nxl) vectors, respectively. Then the equa- 
tions corresponding to u>(t) can be separated from Eq. (10) in the form 


w(t) - I^WfeCO) + [* 2 -(t)]i'(0) + [* 3 (t)Jq/(o) 

+ [B(t-x)][u*(x) + f*(x)]dx 

where [B(t)] is the (3x3) upper left hand corner sub matrix of [4>(t) ] [B^j 

It should be noted that previous investigations [2 ,3] were concerned main- 
ly with the determination of [$, (t) ] and [B(t)] and then with the approximation 
of Eq. (12) by ■ 


( 12 ) 

1 
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u(t) - [^(tJWO) + /* [B(t-T)][u*(T> +l*(T)]dT (13) 

It is now clear that a control system based on Eq. (13) will be successful If 
and [$^(t)] approach zero sufficiently fast. On the other hand, 

as explained in Ref. [1] , such large dancing rates may create fast error growth 
rates in the presence of random disturbances. Hence it is advisable to use 
zero damping rates and obtain the control system from Eq. (12) while measuring 
^ (0) and (0) . 

NUMERICAL SIMULATION 

To illustrate the differences in the responses of the flexible satellite 
under the above mentioned control policies, the results of several numerical 
simulations are now presented. The assumed satellite geometry, mass, inertia, 
and stiffness properties are: 

Nominal spin rate of the satellite m 0.05 rads/sec. 

Nominal principal moments of inertia of the satellite: 

I = 100.0 ; I « 100.0 ; I « 200.0 slug-in, 2 
x 9 y 9 z & 

Total mass of the satellite * 50.0 slugs. 

Diameter of the rigid core « 50.0 ins. 

Length of the beams =30.0 ins. 

Linear mass- density of the beams = 0,0625 slugs/in. 

Area moment of inertia of the beam cross-section ** 0.5 in. 

Young’s modulus for the beam material = 30x10** psi. 

Mass at the ends of the beams * 0.005 slugs. 

Mass in the spring-mass -damper systems = 0,005 slugs. 

Spring constant of the spring-mass- darker systems «= 50.0 lbs/in. 

Damping constant of the s.m.d. systems = 0,1 lbs/sec/In, 

Radial distance of the s.m.d. systems from the center of the rigid core * 
15 . 0 ins . 

Initial values of wCt) : ; 

ui^(0) * 0.01 ; w 2 (0) * 0.0 ; 0)^(0) « 0.01 rads/sec. 

Initial values of (t) , (t) : 

q^(0) -* qj<0> - 0.01 ; qj(0) * 0 , j f 1,2. 
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q^(0) - q£(0) - 0.01 ; q’(0) « 0 , j f 1,2. 

The conplete numerical experiment is performed through the following 
steps; 

Step Is A time interval [0,T] in which the controls are to be effected 
is fixed. In this case T was taken as 5.0 secs. 


Step 2; The satellite is assumed to be rigid, and without controls, such 
that w(t) is given by the solution tr-( t) , of the equation 


[pjji 1 ^) + [p^u 1 ^) + f*(c) - o 


( 14 ) 


Equation (14) is integrated and the responses a)^(t) and a)^(t) are 
plotted in Figure 2. 

Step 3: The satellite is assumed to be rigid and subjected to a time-op- 

timal ’bang-bang’ control, u^(t), such that u^t) is given by the 
solution u*(t) , of the equation 


[P’ ]^ 2 (t) + [P’]aJ 2 (t) + u^t) + f*(t) - 0 


(15) 


1 2 
The ii (t) are computed so as to yield (T) * 0 by the method 

shown in Appendix A. Equation (15) is integrated and the res- 
2 2 

ponses rn^(t) and ai^Ct) are plotted in Figure 3. 

Step 4: The satellite is assumed to be flexible, without control and 

with cl' ( 0) « c^* (0) * 0^, such that u^(t) is given by . Here 


u. 3 (t) - [yt) ]<-)«)) + [B(t-T))f*(T)dT . 


( 16 ) 


3 3 

The responses ^(t) and o^Ct) fromEq. (16) are plotted in Figure 4. 

Step 5: The satellite is assumed to be flexible, with^'(O) — ^'(0) * 0, 

The satellite is subjected to the control torque u^-(t) computed 
in Step 3, such that m(t) is given by uA(t) , where 


Jo 4 (t) «■ t^(t)]w(0) + sl [B(t-T)][f*(x) + u 1 (x)jdT 


( 17 ) 


The responses 0 K(t) and w^(t) from Eq. (17) are plotted in Figure 
5. A 2 


Step 6: The satellite is assumed to be flexible, with * (0) « (0) - jO* 

and subjected to a time-optimal ’bang-bang’ control, u^-(t) , such 
that u(t) is given by oj5(t), where 


Jit) = [yt)]m(0) +/q [B(t-t)]tf*(x) + u 2 (T)3dT 


(18) 
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2 5 

The n (t) are computed so as to yield or (T) * 0^ by the method 

shown in Appendix A. The responses w^(t) and * rom Eq. (18) 

are plotted in Figure 6. 

Step 7: The satellite is assumed to be flexible, with ^ (0) 4 0 1* jf 1 (0) 

and without control, such that &>(t) is given by rn^(t), where 

« 6 (t0 = [4^(0 ]u(0) + [* 2 (t)]i'(0) + [* 3 (t)]i'(0) 

+ /* [B(t-T)][f*(t)]dT . (19) 

The responses (t) end ui^(t) are plotted In Pigure 7. 

Step 8: The satellite is assumed to be flexible, with 3/ (0) 4 0 4 <1/(0) 

and subjected to the control torque u^(t) computed in Step 6* 

such that u^(t) is given by u^(t) , where 

«> 7 (t) - t\(t)iu(o) + t* 2 (t) li'(o) + [* 3 (t)]a.*<o> 

+ /q [B(t-T)][u 2 (T) + £_*(t) ]dt . (20) 

7 7 

The responses W^(t) and w^Ct) are plotted in Figure 8, 

Step 9: The satellite is assumed to be flexible, with 3/ (0) 4 0 4 fi* (0) . 

Xt is also subjected to a time-optimal 1 bang-bang 1 control, 

3 . r 

u (t) , such that w{t) is given by uT'(t) y where 

uV) - {* x (t) ]u(o) + (* 2 (t)]i'(o) + [* 3 (t)ia.’(o) 

+ /* [B(t-T)Uu 3 (T) + £*(T)]dT . (21) 

3 8 

The torques _u (t) are also computed to yield id (T) = by the 

8 9 

method shown in Appendix A. The responses o]^(t) and a>^(t) from 
Eq. (21) are plotted in Figure 9. 

COMPARISON AND EVALUATION 

One important result of the simulation, as seen from Figures 2 and 3, is 
that the control sequence uvt) is very effective on the rigid model of the sa- 
tellite. But Figure 5 shows that, for the same values of m(0), u^(t) produces 
unwanted non-zero values of tjj(T) when it is applied to the flexible satellite mo- 
del, although c^(Q) and ^(0) are assumed to be zero. Thus, another important 

0 

result, presented in Figures 5 and 6, shows that ju (t) is more effective than 
u^t) when a flexible satellite model is considered. Up to this point, then, we 
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have essentially the same conclusion as that in Kefs. [2], [3], that for a flex- 
ible satellite the control should not be based on a rigid model. The difference 
between Refs. [2], [3] and the present study is in the adopted control policy, 
'Bang-bang' controls have been used here instead of linear proportional control. 

The most inportant results are presented in Figures 8 and 9. When the 
^(0) and ^(0) are observed and found different from zero, u^(t) does not lead 
to the required *ero values of m(T) . In contrast, u3(t), which is based on the 
observed values of ^(0) and <^(0) , yields zero values of jd(T) . Another point 
to be considered is the divergence of ui(t) from zero in the two cases. The 

maximum divergence of jtn(t) and u^(t) is 11.0 x 10 ^ rads/sec, while that with 
3 -4 

ju (t) is only 7.0 x 10 rads/sec. This bears out the theoretical claims that 
a control based on Eq. (12) is more effective than one based on Eq , (13) and 
that the effectiveness of a control system can be greatly Improved if the de- 
flections of the flexible elements of a satellite are observed. 
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APPENDIX A 



The method of computing the time-optimal control torques for a system 
given by 

x(t) » 2 _(t) + /*!K(t-x)]u(x)dx (22) 

is now presented. Reference [5] presents computing algorithms and other de- 
tails of the method. In Eq. (22), x(t) is the output vector of the system, 
u(X) is the control vector, and ^(t) and [K(t)J are known vector and matrix 
functions of the time, t. 

It is assumed that, for a given t = T, u(t) should be such that 

_x(T) > 0 (23) 

and |ju(t) | for all t is a minimura. Thus, the minimum time problem is converted 
to the equivalent minimum control effort problem. The solution for u(t) is 
then given by [5]. 

UjXt) * U(T)y* sgn [I AjK^T-t)] (24) 

where 

U(T) = 1.0/ [min £ |Z A*K . (T-x) | dx] (25) 

A* i U i 1 13 

such that 

S A*y i ( T) = 1.0 (26) 

and 

P* - - \ [KjTiCr) . (27) 

[K l ] lj “ A [K lj (t_T)] sgn ^*K rj (T-x)]dx (28) 

The summation convention of repeated indices is not to be used in Eqs v (24) 
to (27) above. 


79 



• SPRING-MASS 
-DAMPER SYSTEM 

^ — RIGID BODY 

— 4 — — — 


CO-ORDINATE 

AXES 


BEAM NO* 3— 1 
END -MASS - 

BEAM NO. 2 


SPRING- MASS 
-DAMPER SYSTEM 


BEAM NO. I 


BEAM NO. 3 





1 

BEAM NO. 4 

FIGURE 1 

0 

: SATELLITE CONFIGURATION 








(_) 

LU 

CO 

\ 


cz 

or 


to 

UJ 


o 

_J 

LU 

> 


CO 


* 10 * 


RIGID RESPONSE 
WITHOUT CONTROL 


MC-* 



10-0 
—I 7.5 
5.0 
2.5 
.0 

-1 -2 -S 

-5.0 
-7-5 
10.0 


* 10 -* 


• 00 .50 1.00 1.50 2. CO 2-50 3-00 3.50 

TIME , EEC. 


4.00 4.50 S.05 


- 1*0 
1C* 4 


F igure 2 


RNG. VELOCITIES . CRD. /SEC- 


! 

I 

| 


RIGID RESPONSE 
WITH CONTROL 


•00 .50 1-00 1-50 2-00 2-50 3-00 3-50 4.00 4.50 5-00 



TIME 


C" r~ 


Figure 3 


'PNG. VELOCITIES , RRO./SEC- 




PNG. VELOCITIES . RRD./SEC. 


FLEXIBLE RESPONSE 
WI l H RIGID CONTROLS 



TIME 


Figure 5 


•RNG . VELOCITIES , RflO./SEC. 


F L E X I B'L E RESPONSE 
FLEXIBLE CONTROL 


.00 -50 I .00 1.50 2-00 2-50 3-00 3-50 4-00 4-50 5-00 



TIME , SEC 

Figure 6 


■RNC. VELOCITIES . K’HD./SEC- 


OBSERVED' RESPONSE 
WITHOUT. CONTROL 

J 

.00 -50 1 -00 1-50 Z-00 Z -50 3-00 3-50 4-00 4-50 5-00 



TIME , SEC- 


Figure 7 










PNG. VELOCITIES , RHD./SEC. 


OBSERVED RESPONSE 
OBSERVED CONTROL 


.00 -50 1 -00 1-50 2-00 2-50 3-00 3-3C 4 - CO ^-50 5-00 



TIME , EEC * 

Figure 9 


• 0 

.< 

0 

0 

0 

,0 

,0 

.0 

0 

■ 0 

,0 

.0 


